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For  the  past  five  years  the  Buoy  Group  at  Woods  Hole  Oceanographic 
Institution  has  included  temperature  as  one  of  the  variables  recorded 
in  its  current  meters.  These  measurements  began  with  the  first  success- 
ful deployments  of  Vector  Averaging  Current  Meters  (VACMs)  in  1971. 
Circuitry  designed  for  making  highly  accurate  temperature  measurements 
has  been  included  in  all  the  Buoy  Project's  VACMs.  During  the  past  year 
we  have  begun  to  add  similar  circuitry  to  the  EG&G  850  current  meters. 

This  report  is  intended  to  describewhat  we  have  learned  about  making 
water  temperature  measurements  with  VACMs. 

Among  the  authors,  K.  Schleicher  and  A.  Bradshaw  are  responsible 
for  thermistor  calibrations,  J.  Dean  for  quality  control  of  VACM  main- 
tenance and  calibration  work,  and  R.  Payne  for  data  analysis. 
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1 . Introduction 

The  Vector  Averaging  Current  Meter  (VACM)  was  designed  several  years 
ago  by  engineers  at  the  Woods  Hole  Oceanographic  Institution  in  response 
to  the  needs  of  scientists  associated  with  the  Buoy  Project  at  that 
Institution.  After  the  first  few  prototypes  were  built  and  tested  the 
design  was  turned  over  to  the  AMF  Electrical  Products  Division  of  Alexandria, 
Virginia  which  builds  and  sells  the  instrument.  Included  in  the  original 
design  and  in  the  standard  instrument  sold  by  AMF  is  circuitry  capable  of 
making  accurate  temperature  measurements.  This  report  describes  calibration 
techniques  we  have  found  effective  for  both  sensor  and  electronics  as 
well  as  evaluations  of  the  performance  of  the  two  components  independently 
and  as  a system. 

System  description 

A thermistor  is  mounted  in  the  lower  plate  of  a chassis  assembly 
which,  in  turn,  is  mounted  on  the  end  plate  of  the  pressure  housing  of 
the  VACM.  An  electrical  current  through  the  thermistor  is  sensed  by  an 
electronic  circuit  which  provides  a series  of  pulses  whose  rate  is  pro- 
portional to  the  current.  These  pulses  are  summed,  typically,  over  a 
15  minute  interval  and  the  total  recorded  as  part  of  the  VACM  data  record 
on  a magnetic  tape  cassette.  This  arrangement  allows  true  time  averaging 
of  the  thermistor  resistance  over  the  entire  recording  interval.  For 
most  oceanographic  applications  this  is  equivalent  to  a time  average  of 
temperature.  This  technique  has  been  used  since  about  1971  and  has 
proved  to  be  quite  reliable  and  accurate. 

The  standard  VACM  uses  a Yellow  Springs  Instrument  Corp.  (YSI) 
thermistor,  part  #44032,  with  0.1°  C interchangeability.  In  addition  to 
these  we  have  recently  been  using  in  our  routine  operations  some  Ther- 
mometries Corp.  thermistors  with  quite  similar  characteristics  developed 
for  the  Internal  Wave  Experiment  (IWEX) . This  report  will  describe  results 
on  the  YSI  thermistors  only  since  we  have  accumulated  the  most  experience 
with  them.  The  epoxy  encapsulated  thermistors  are  potted  in  an  aluminum 
screw  for  attachment  to  the  VACM  chassis.  The  nominal  resistance  of  the 
thermistor  is  30,000  ohms  at  25°  C with  a temperature  coefficient  of 
-4.5  %/°  C. 
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The  thermal  time  constant  (time  required  for  the  thermistor  to 
sense  63%  of  a step  change  in  temperature)  for  the  device  when  installed 
in  a VACM  end  cap  is  a function  of  water  current  speed  and  has  been 
measured  to  be  approximately  100  seconds  for  an  average  current  of 
14  cm/sec. 

We  believe  that  the  absolute  accuracy  of  our  temperature  measurements 
in  VACMs  is  better  than  .01°  C.  This  has  been  substantiated  by  the 
intercomparisons  described  in  this  report  and  by  analysis  of  the  errors 
inherent  in  the  the  instrument  (Section  6) . Thermistor  stability  and 
the  practical  difficulties  of  making  extremely  precise  temperature  calib- 
rations appear  to  be  the  limiting  factors.  Making  a series  of  calibrations 
over  several  years  and  culling  unstable  thermistors  are  necessary  for 
highest  accuracy. 

We  might  point  out  that  making  meaningful  absolute  temperature 
measurements  with  accuracies  of  .01°  C requires  stability  and  accuracy 
of  depth  which  strains  present  day  mooring  techniques.  In  the  main 
thermocline  in  the  Atlantic  Ocean  the  temperature  gradient  is  of  the  order 
of  20  m°  C/m.  A temperature  accuracy  of  10  m°  C there  is  equivalent  to 
a depth  accuracy  of  50  cm. 

There  are  other  reasons  for  making  precise  calibration  of  sensors 
and  associated  circuitry.  Currently  many  of  our  moorings  are  out  for 
periods  approaching  1 year.  Much  of  the  use  made  of  temperatures,  such 
as  computing  low  frequency  heat  fluxes,  depends  more  on  the  stability 
of  temperature  measurements  than  on  their  absolute  accuracy.  A drift  of 
20  m°  C/year  in  a thermistor  can  make  temperature  measurements  well  below 
the  main  thermocline  virtually  meaningless  unless  it  can  be  corrected  for. 
Accuracy  of  drift  rate  estimates  depends  on  the  uncertainty  in  individual 
calibrations  and  the  length  of  time  over  which  the  calibrations  have  been 
made.  Doubling  the  uncertainty  in  accuracy  of  individual  thermistor 
calibrations  doubles  the  length  of  time  over  which  calibrations  must  be 
made  to  achieve  the  same  accuracy  in  drift  rate  estimates.  We  find  that 
a minimum  of  3 calibrations  over  2 years  are  required  to  determine  drift 
rates  to  our  specifications.  Although  highly  accurate  calibrations  are 
expensive,  an  increase  in  calibration  period  would  leave  us  with  virtually 
no  certifiably  stable  thermistors. 
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2.  V/'F  Converter  Operation 

In  the  VACM  the  thermistor  is  used  to  vary  the  electrical  current 
into  a converter  (called  the  voltage  to  frequency  (V/F)  converter). 

Since  the  period  of  the  pulse  train  output  of  the  converter  is  propor- 
tional to  the  input  current,  the  average  resistance  (and  therefore  tem- 
perature) of  the  sensor  over  a precisely  measured  time  interval  can  be 
obtained  by  counting  the  pulses  in  that  interval.  The  pulses  are  summed, 
typically  for  15  minutes.  A technical  discussion  of  the  circuitry  can 
be  found  in  the  VACM  technical  manual,  AMF  Publication  #SLS  106-11419, 
Section  2.2.1.  The  technique  described  has  been  used  since  about  1971 
and  has  proved  to  be  quite  precise  and  reliable. 

Referring  to  Figure  1,  operational  amplifier  (op  amp) , A,  and 
feedback  capacitor,  C , form  a current  integrator.  Ideally,  an  op  amp 
has  infinite  gain,  the  input  summing  points  (+  and  -)  have  zero  potential 
between  them,  and  no  current  can  flow  into  the  input  terminals.  A stable, 
precise  reference  voltage,  V,  is  applied  to  Rt#  the  series  combination 

of  thermistor  R and  a fixed  precision  resistor  R , resulting  in  a 
t s 

current  i = V/R  . Since  no  current  can  flow  into  the  op  amp,  current  i 
r T r 

must  flow  in  C^  and  equal  if. 


Figure  1 
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V/F  Converter 
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Periodically,  a precise  electrical  charge,  Q,  is  applied  to  the 
op  amp  input  summing  point  and  is  transferred  to  capacitor  . This 
charge  Q is  dispensed  by  if,  and  the  time  period,  P,  required  to 
remove  the  charge  is  a function  of  the  current  i^. 

The  charge  Q accumulated  on  the  precision  capacitor  C by  reference 
voltage  V is 


Q = CV. 

K 

. 

The  input  current  flowing  in  the  thermistor  is 


ir  - v/V 

The  discharge  of  Q by  current  i^  in  time  P is 
if  = Q/P  = CV/P. 

Since 


then 


V/R^  = CV/P,  or 

P = R^C  = (Rt  + Rs)C. 

The  period  P is  proportional  to  the  total  resistance  Rj,,  hence 
related  to  the  thermistor  resistance  Rfc.  Rfc  is  converted  to  temperature 
by  substitution  into  the  calibration  equation  (see  Section  5) . 

Solid  state  FET  switches  transfer  the  charge  from  C to  via 

op  amp  A.  After  an  appropriate  interval  (nominally  400  microseconds)  the 
switches  reconnect  C to  V to  recharge  C.  continues  to  discharge 

and  after  time  interval  P,  is  discharged,  the  switches  again  trans- 

fer the  charge  on  C to  the  integrator,  and  the  cycle  begins  again.  C 
and  C2  form  a voltage  divider  to  reduce  the  voltage  at  the  input  of  A 
to  prevent  driving  the  op  amp  input  into  saturation  where  it  would  draw 
current. 
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Additional  circuitry  controls  the  switching  and  provides  closed 
loop  operation  whereby  the  cycle  is  reinitialized  as  soon  as  the  charge 
is  dispensed.  The  circuit  also  provides  a slight  delay  in  the  switching 
(about  10  micro-seconds)  to  allow  break-before-make  action  of  the  switches 
to  prevent  loss  of  charge. 

Since  the  frequency  output  is  a function  of  the  total  resistance  Rt 
and  the  capacitance  of  capacitor  C (P  = R^C) , the  circuit  is  relatively 
insensitive  to  variations  in  components  in  the  circuit  except  for  Rg 
and  C.  The  series  resistor  Rg  is  chosen  to  produce  a nearly  linear 
relationship  between  the  current  in  the  thermistor  and  its  temperature. 
With  this  resistor  value  equal  to  37,300  ohms,  the  current  in  the  ther- 
mistor-resistor network  is  linear  with  temperature  within  1.7%  over  the 
range  of  0°  C to  30°  C.  Although  the  network  is  somewhat  non-linear 
temperature  excursions  in  the  ocean  do  not  often  cause  biases  in  a 
15  minute  average.  The  output  frequency  of  the  converter  is  proportional 
to  current  input  within  .01%. 

There  is  a slight  error  in  the  period  (not  exactly  equal  to  R^C) 
and  is  probably  due  to  a finite  reaction  time  of  the  op  amp  to  the  charge 
applied  at  its  input.  An  empirical  constant,  K,  can  be  determined  from 
tests,  resulting  in  a more  nearly  exact  equation 


P = Rt  C + K. 

K has  been  found  to  be  about  0.5  microseconds  for  a typical  instrument. 

The  expression  for  computing  thermistor  resistance  from  the  period  of 

the  V/F  converter  output  signal  is  straightforward  to  derive,  if  rather 

involved.  One  begins  with  the  valid  assumptions  that  the  voltage  across 

Rfc  and  R^in  Fig.  1 is  constant  and  that  the  frequency  of  the  V/F  converter 

is  directly  proportional  to  the  current  through  R and  R . The  resulting 

t s 

expression  is: 


,1  x 

F + K2  s 

Ki  - K2  - 1 
P 


where  P = T^/N  = period  of  V/F  output 


Ki  ■ <i  - 5 l/D 

2 1 


K2  " 


P2(RS+V  P1(RS+R2) 


/D 


D = 


RS  + R2  RS  + R1 


and  R^,  R^»  P^»  are  t^le  input  resistances  and  output  periods  of  two 
V/F  calibration  points. 
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3.  V/F  Converter  calibration 

Since  the  V/F  converter  transforms  the  resistance  changes  of  the 
thermistor  into  electrical  signals  it  is  essential  that  it  be  calibrated 
as  carefully  as  the  thermistor.  Our  calibration  procedure  is  designed 
to  monitor  the  long  term  behavior  as  well  as  provide  a calibration  for  a 
specific  setting  of  an  instrument. 

Before  and  after  each  deployment  the  V/F  converter  is  calibrated  at 
seven  precise  (.001%)  input  resistances  corresponding  to  nominal  thermistor 
resistance  values  at  seven  temperatures.  The  values  are  shown  in  Table  1. 
The  period  at  the  5°  C point  is  adjusted  to  read  1820.4  ys  at  74434 
and  the  seven  periods  are  recorded  again,  both  at  room  temperature  and  in 
a cold  room  at  V5°  c.  A record  is  made  at  the  periods  before  and  after 
adjustment  which  allows  us  to  monitor  the  temporal  stability  of  each 
( r'  converter.  The  cold  room  values  are  used  in  decoding  most  data  since 

they  correspond  closely  to  the  in  situ  conditions  of  the  instruments  in 
typical  moorings. 

I 

I Table  1 


V/F  Converter  Calibration  Results 


Temperature 

Input  Resistance 

Period  (x  10 

-6 

seconds 

0°  C 

04972  fi 

2154.871 

+ 

.075 

5 

74434 

1820.412 

+ 

.004 

10 

58747 

1564.911 

+ 

.056 

15 

46675 

1368.2 

20 

37299 

1215.640 

+ 

.093 

25 

29998 

1096.765 

+ 

.109 

30 

24269 

1003.464 

± 

.127 

The  period  values  in  Table  1 are  the  means  and  standard  deviations 
of  the  data  points  at  each  temperature  for  the  175  calibrations  investigated. 
Only  a nominal  value  has  been  included  for  the  15°  C point  because  we 
discovered  some  of  the  calibrations  were  made  with  an  incorrectly  labeled 
precision  resistor. 


The  period-resistance  relationship  is  linear.  In  decoding  the  data 
from  a current  meter  we  use  the  5°  C and  25°  C points  to  determine  the 
relation.  We  checked  the  linearity  for  175  calibrations  of  45  instruments 
by  comparing  the  actual  calibration  points  and  the  values  computed  from 
the  straight  line  determined  by  the  5°  C and  25°  C points.  The  result 
was  that,  in  all  but  five  calibrations,  the  error  introduced  by  using  the 
straight  line  was  equivalent  to  1 m°  C and  in  those  five  it  was  equi- 
valent to  about  2 m°  C inaccuracy  in  the  decoded  temperatures. 

Using  the  5°  C and  25°  C points  from  Table  1 the  period-resistance 
relation  is 


PER  (ys)  = 608.243  + 0. 01628515*R. 

A least  squares  fit  of  all  but  the  15°  C point  yields 

PER  (ys)  = 608.231  + 0. 01628516*R. 

The  two  equations  fit  the  data  equally  well.  The  standard  deviation  of 
the  first  is  0.013  ysec,  and  of  the  second,  0.018  ysec.  The  accuracy  in 
period  required  for  ±1  m°  C accuracy  in  temperature  is  ±0.08  ysec  at  0°  C, 
±0.05  ysec  at  10°  C,  and  ±0.03  ysec  at  20°  C. 

We  have  also  fit  an  equation  like  the  first  of  the  preceding  para- 
graph to  the  data  of  the  175  individual  V/F  calibrations.  After  sol- 
ving each  equation  for  R and  subtracting  the  input  resistances  we  com- 
puted an  r.m.s.  residual  for  each  calibration.  The  mean  of  this  residual 
over  all  the  calibrations  is  0.9  ohms.  All  but  three  of  the  residuals 
were  less  than  3 ohms.  One  ohm  represents  0.22  m°  C at  0°  C and  0.99  m°  C 
at  30°  C for  our  standard  thermistors. 

The  calibrations  also  show  that  the  stability  of  the  V/F  con- 
verters is  such  that  they  will  maintain  1 m°  C accuracy  over  the  period 
of  one  year,  a stability  of  3 parts  in  105. 
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4.  Thermistor  Calibrations:  Methods  and  Facilities 

All  of  our  thermistor  calibrations  have  been  made  with  the  same 
variety  of  equipment:  a constant  temperature  bath,  a platinunr  resistance 

thermometer  (PRT)  and  bridges  to  measure  the  resistance  of  both  the  PRT 
and  the  thermistor.  Calibrations  began  in  January  1973.  In  August  1973 
(just  prior  to  batch  11)  there  was  a bath  modification  resulting  in  an 
improvement  in  accuracy.  In  December  1974  (just  prior  to  Batch  45)  most 
of  the  system  was  replaced  for  some  improvement  in  accuracy  and  a large 
improvement  in  convenience  a no  speed . 

We  will  describe  both  calibration  systems  and  the  errors  connected 
with  them. 

A.  System  1.  January  1973  - December  1974.  Batches  1-10 

Figure  2 shows  a diagram  of  the  first  bath.  It  was  one  used  pre- 
viously (Bradshaw  and  Schleicher,  1970)  but  its  precision  of  control 
was  improved  for  this  application.  It  consisted  of  an  outer  and  an 
inner  bath,  both  filled  with  a stirred  water-ethylene  glycol  solution. 

The  outer  bath  had  heating  and  cooling  coils.  The  inner  bath  was  heated 
to  maintain  it  1-2°  C above  the  temperature  of  the  outer  bath.  Within  the 
inner  bath  was  a kerosene-filled  aluminum  vessel  containing  the  thermistor 
mounting  block  (see  Fig.  3).  The  kerosene  was  stirred. 

The  thermistor  mounting  block  was  made  from  a 10"  length  of  aluminum 
hexagonal  bar  stock  1 1/4"  across  the  flats.  The  thermistors  were 
screwed  into  radial  tapped  holes  in  the  flat  sides,  4 to  a side.  The 
platinum  resistance  thermometer  was  mounted  in  a 3/8"  diameter  hole  drilled 
to  within  1/2"  of  the  bottom.  Mercury  was  used  in  the  bottom  of  the  well 
for  good  thermal  contact. 

A guarded  Wheatstone  bridge,  Leeds  and  Northrup  type  4737-A20, 
was  used  to  measure  thermistor  resistance. 

A Leeds  and  Northrup  type  8163-C  platinum  resistance  thermometer 
was  used  to  measure  the  temperature  of  the  thermistor  mounting  block. 

The  resistance  of  the  thermometer  was  measured  with  a Leeds  and  Northrup 
type  8068-B  Mueller  bridge  using  a nanovoltmeter  as  a null  meter.  The 
calibration  of  the  PRT  was  checked  periodically  at  the  triple  point 
temperature  of  water  with  a Trans-Sonic  Equiphase  type  130  triple  point 
cell . 
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Up  to  23  thermistors  were  mounted  in  the  aluminum  block.  The 
thermistor  leads  were  connected  to  cables  going  to  the  Wheatstone  bridge 
via  terminal  strips.  The  24th  terminal  pair  of  these  strips  was  usually 
connected  to  a pair  of  bare  wires  immersed  in  the  kerosene.  The  resis- 
tance of  this  open  line  was  monitored  with  a megohmmeter  as  a check  on 
any  condition  (contaminated  oil,  condensation,  etc.)  which  could  cause 
a significant  shunting  error  in  the  thermistor  resistance  measurements. 

Resistance  of  the  thermistors  was  measured  in  sequence  at  a set  of 
temperatures  starting  at  approximately  06  C and  increasing  in  5°  C steps 
to  30°  C.  Voltage  across  the  thermistor  was  set  to  correspond  to  that 
appearing  in  the  current  meter  to  avoid  differences  in  power  dissipated 
in  the  thermistor.  The  whole  sequence  of  measurements  took  2 or  3 days. 
Estimated  Random  and  Systematic  Errors  of  Measurements 

During  this  period  (January  1973  - January  1975)  a change  was  made 
in  the  apparatus  which  affected  the  estimated  precision  of  measurement. 
This  occurred  in  August  1973  when  a temperature  gradient  was  discovered 
in  the  oil-filled  vessel.  From  tests  made  at  this  time  (by  interchanging 
the  thermistors  in  the  highest  and  lowest  position  in  a similar  mounting 
block)  the  largest  temperature  difference  from  the  top  of  the  mounting 
block  to  the  bottom  appeared  to  be  .004°  C at  a bath  temperature  of  0°  C. 
After  changing  the  oil  level  and  introducing  baffles  inside  the  aluminum 
vessel,  retesting  showed  that  the  difference,  if  any,  was  less  than 
.001°  C.  The  estimated  random  error  prior  to  Batch  11  (September  7, 

1973)  includes  a factor  based  on  the  larger  gradient. 

The  temperature  variation  in  the  water  bath  due  to  the  control 
cycle  (-  1 cycle/min)  was  usually  less  than  1 m°  C.  The  oil-filled 
vessel  further  reduced  the  cycling  variation  to  the  order  of  several 
tenths  of  a millidegree  as  measured  by  the  platinum  thermometer.  The 
longer  term  temperature  variation  inside  the  mounting  block  could  be 
followed  by  the  thermistors  and  by  the  platinum  thermometer  and  under 
typical  conditions  did  not  vary  by  more  than  1 m°  C over  the  time  required 
to  make  all  measurements  at  one  temperature  point  (1  to  1 1/2  hour  for 
23  thermistors) . 


12 


The  values  of  the  "a"  and  "6"  constants  of  the  platinum  ther- 
mometer are  traceable  to  a NBS  calibration  (see  NBS  Monograph  126,  1973, 
for  symbol  nomenclature)  and  its  resistance  was  checked  regularly  at  the 
triple  point  of  water.  This  resistance  changed  by  the  equivalent  of 
about  2.5  m°  C over  the  two  year  period.  The  Mueller  bridge  calibration 
corrections  are  traceable  to  a NBS  calibration.  This  thermometer-bridge 
combination  was  checked  in  February  1975  against  that  of  a newer  Leeds 
and  Northrup  platinum  thermometer  and  Guildline  Current  Comparator 
Resistance  Bridge  (see  section  B in  this  chapter) . The  temperatures 
indicated  by  the  two  systems  differed  by  not  more  than  0.7  m°  C at  30°  C. 

The  warranted  (uncorrected)  accuracy  of  the  type  4737  guarded 
Wheatstone  bridge  used  to  measure  the  thermistor  resistances  is  ±100  ppm 
at  25°  C (equivalent  to  ±2  m°  C) . The  accuracy  is  derated  by  10  ppm  per  °C 
difference  from  25°  C.  Corrections  to  the  bridge  dial  readings  were 
found  at  25°  C using  resistance  ladders,  one  element  of  which  had  been 
measured  to  an  accuracy  of  5 ppm.  The  accuracy  of  the  corrected  guarded 
Wheatstone  bridge  readings  at  25°  C was  then  taken  as  about  5 ppm,  equi- 
valent to  0.1  m°  C when  applied  to  the  measurement  of  thermistor  resistance 

The  corrections  described  above  were  used  throughout  the  January 
1973  to  December  1974  period.  In  November  1974,  they  were  checked  by 
measuring  the  resistances  of  the  Guarded  Wheatstone  Bridge  using  a 
Guildline  Comparator  Bridge  and  Guildline  standard  resistors.  Differences 
of  up  to  20  ppm  (0.4  m°  C)  were  found.  These  could  be  accounted  for  by 
the  temperature  of  the  Guarded  Bridge  (23-24°  C instead  of  25°  C)  and  by 
a drift  in  the  Guildline  standard  resistors;  however,  this  figure,  0.4  m°  C 
instead  of  0.1  m°  C as  found  earlier,  is  taken  as  the  systematic  error 
associated  with  the  type  4737  bridge. 

The  first  two  columns  of  Table  2 sum  up  our  estimates  of  the  errors 
inherent  in  the  VACM  thermistor  calibration  for  the  calibration  system  in 
use  until  January  1975. 

B.  System  II.  January  1975  (Batch  45)  to  present 

In  January  1975  we  began  using  a new  set  of  thermistor  calibration 
apparatus  which  has  moderately  improved  the  calibration  accuracy  and 
substantially  increased  the  speed  and  convenience  of  making  calibrations. 
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The  new  bath  system  consists  of  a Tronac  Model  400  constant  tem- 
perature bath  with  a Tronac  Model  40  temperature  controller.  This  is 
a single  bath  system  making  it  easier  to  change  the  bath  temperature 
rapidly. 

The  new  thermistor  mounting  block  (Fig.  4)  was  machined  from  an 
aluminum  cylinder,  4 inch  diameter  by  4 inch  long.  Six  lengthwise, 

1 inch  deep  cuts  around  the  cylinder  allow  recessed  mounting  of  4 thermistors 
per  cut.  The  platinum  thermometer  is  mounted  in  a 3/8  inch  diameter  hole 
drilled  through  the  axis  of  the  cylinder.  The  mounting  block  is  enclosed 
in  an  aluminum  housing  which  is  purged  with  dry  nitrogen  and  hermetically 
sealed  before  being  placed  in  the  bath. 

The  new  platinum  resistance  thermometer  is  a Leeds  and  Northrup 
Model  8167-25B.  Its  resistance  is  measured  with  a Guildline  Instruments 
Model  9975  thermometer  bridge.  As  used  ’ n thermometry,  the  bridge  de- 
termines the  ratio  of  the  platinum  probe  resistance  and  a high  precision, 
stable  resistor  mounted  in  a temperature  controlled  air  bath.  As  with 
the  Mueller  bridge,  a four  terminal  measurement  of  the  probe  resistance 
is  made  and  the  effect  of  the  thermal  emf's  is  eliminated  by  current 
reversals.  Readings  can  be  made  in  far  less  time  and  with  greater  precision 
than  with  the  Mueller  bridge  used  previously. 

The  rest  of  the  equipment,  the  Wheatstone  bridge  and  null  meter  for 
measuring  thermistor  resistance,  and  the  triple  point  cell,  remains  the 
same  as  before. 

The  calibration  procedure  is  the  same  as  before  except  that  the 
amount  of  time  required  to  calibrate  a normal  batch  of  24  thermistors 
is  now  12  hours  or  less. 

Estimated  Random  and  Systematic  Errors  of  Measurements 

During  this  period  (January  1975  to  the  present)  the  random  and 
systematic  errors  of  measurements  were  estimated  to  be  somewhat  less 
than  in  the  previous  period.  They  are  summarized  in  column  3 of  Table  2. 

The  mounting  block  gradient,  one  of  the  main  sources  of  error 
in  precision  during  the  previous  period,  was  improved  through  the  use 
of  the  newer  apparatus.  A test  done  in  March  1976  which  included  a 
careful  measurement  of  the  vertical  profile  of  the  temperature  along 
the  axis  of  the  platinum  thermometer  well  as  well  as  the  exchange  of 


thermistor  position  indicated  a maximum  variation  of  less  than  0.6  m°  C. 

As  a result  of  the  improved  temperature  control  of  the  new  bath 
the  error  due  to  the  short  term  temperature  cycling  of  the  mounting 
block  was  estimated  to  be  not  more  than’O.l  m°  C.  Actually,  much  of  the 
time,  the  temperature  read  by  the  platinum  probe  did  not  change  during 
the  scanning  of  24  thermistors  at  one  temperature  (1/2  to  3/4  hours) 
by  more  than  0.1  or  0.2  m°  C.  When  the  drift  was  larger  than  this  the 
tracking  of  the  recorded  thermistor  temperature  and  platinum  thermometer 
temperature  was  close  enough  so  that  the  errors  were  still  kept  within 
the  0.1  m°  C figure. 

The  largest  source  of  error  in  precision  with  this  n<w  equipment 
was  due  to  the  drift  of  the  platinum  thermometer  between  triple  point 
readings.  Of  course,  most  of  this  error  can  be  corrected  by  interpolating 
linearly  between  check  points  but  since  this  would  have  delayed  the 
processing  and  reporting  of  the  data  it  was  not  done.  In  the  future  it 
is  planned  to  take  more  frequent  triple  point  checks  to  reduce  this 
source  of  error. 

The  random  errors  due  to  temperature  effects  on  the  bridge  used 
to  measure  the  thermistor  resistances  remain  the  same  as  in  the  previous 
period.  There  are  no  significant  temperature  effects  (<  0.1  m°  C)  on 
the  Guildline  bridge. 

The  systematic  error  of  the  platinum  resistance  thermometer  is 
estimated  to  be  0.5  m"  C from  information  obtained  from  the  manufacturer 
Leeds  and  Northrup  (L&N) . 

The  Guildline  bridge  error  in  ratio  measurement  is  certified  by  the 
manufacturer  not  to  exceed  ±(2  parts  in  10^  of  reading  + 1 step  of  last 
dial)  . This  would  cause  an  insignificant  sytematic  error  in  the  tem- 
perature (<  .1  m°  C) . 

The  systematic  errors  associated  with  the  type  4737  bridge  and 
the  triple  point  cell  remain  the  same  as  before. 


Random  and  Systematic  Errors  of  Measurements  in  m°  C Equivalents 
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5.  Thermistor  Calibration  Results 

In  this  section  we  will  describe  what  we  have  learned  from  calibrations 
to  date  of  our  YSI  type  44032  thermistors. 

We  now  own  163  YSI  type  44032  thermistors  which  have  each  had  at 
least  one  calibration.  Of  these,  69  meet  the  criteria  for  being  included 
in  this  study,  i.e.,  at  least  three  calibrations  over  at  least  two  years. 

We  intend  to  continue  regular  calibrations  of  all  our  thermistors. 

Our  analysis  of  thermistor  performance  is  based  on  the  outputs  of 
the  computer  programs  CALIBB  and  CALCOLIN,  described  in  Appendices  II 
and  III.  CALIBB  fits  the  resistance,  temperature  calibration  data  to 
the  expression  (Bennett,  1972,  Steinhart  and  Hart,  1968) 

T_1  = A + BfcnR  + cUnR)3 

by  a least  squares  technique.  A residual  is  computed  for  each  data 
point  as  the  difference  between  the  observed  temperature  and  a temperature 
computed  from  the  observed  thermistor  resistance  and  the  A,  B,  C 
constants  from  that  calibration.  The  r.m.s.  value  of  this  residual  for 
the  seven  data  points  in  a calibration  is  printed  out  in  CALIBB  and  falls 
in  the  range  0-1  m°  C for  nearly  all  of  our  calibrations.  The  mean 
value  for  all  our  past  calibrations  is  .2  m°  C,  which  shows  how  well  the 
above  equation  fits  the  data. 

Using  CALCOLIN  we  can  estimate  the  rate  at  which  the  temperature 
indicated  by  a thermistor  is  changing  relative  to  the  true  temperature  of 
its  environment.  We  should  remember,  however,  that  what  is  actually 
changing  is  the  resistance  of  the  thermistor  at  a given  temperature. 

Table  3 and  Figure  5 sum  up  the  calibration  histories  of  the  69 
thermistors  with  sufficient  calibration  history.  Table  3 lists  for  each 
thermistor,  in  order  of  increasing  drift  rate  magnitude,  the  total  number 
of  calibrations,  the  elapsed  time  in  years  between  the  first  and  latest 
calibration,  the  mean  drift  rate,  and  the  standard  deviation  about  the 
mean  drift.  The  drift  rate  and  standard  deviation  are  means  over  the 
least  squares  fit  at  each  of  the  seven  calibration  temperatures  and  are 
computed  in  CALCOLIN.  In  this  table  drift  rate  has  the  opposite  sign 
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from  the  number  in  the  CALCOLIN  output  since  it  is  the  rate  at  which 
temperature  would  appear  to  change  if  the  thermistor  were  kept  immersed 
in  a bath  whose  temperature  was  held  constant.  Positive  drift  rate 
indicates  that  thermistor  resistance  at  constant  temperature  is  decreas- 
ing with  time. 

The  standard  deviations  listed  in  Table  3 show  that  our  calib- 
rations are  consistent  with  each  other  to  better  than  1 m°  C.  Standard 
deviations  are  less  than  this,  for  the  most  part,  except  for  the  higher 
drift  rates.  There  is  other  evidence  for  the  repeatibility  of  our 
calibrations.  Calibration  batches  52  and  53  were  done  within  1 week 
of  each  other  and  19  thermistors  were  common  to  both  batches.  The  mean 
difference  between  the  two  calibrations  for  each  of  the  19  thermistors 
was  .3  m°  C.  The  r.m.s.  difference  was  .97  m°  C. 

Figure  5 is  a histogram  of  the  mean  drift  rate  magnitudes  for  the 
69  thermistors  in  Table  3.  The  dividing  line  we  have  chosen  between 
stable  and  unstable  thermistors  is  arbitrary  but  can  be  rationalized. 

The  absolute  accuracy  of  our  temperature  measurements  is  no  worse  than 
±10  m°  C and  may  be  as  good  as  ±5  m°  C.  For  a typical  deployment  of 
9-12  months  we  would  prefer  that  the  thermistor  drift  be  less  than 
5 m°  C so  that  correction  for  drift  is  not  necessary.  The  dip  in  the 
distribution  around  3 m°  C/year  provides  a natural  place  to  draw  the  line. 

Since  the  thermistors  we  use  are  standard  off-the-shelf  units 
with  no  preselection  applied  we  would  expect  the  distribution  shown  in 
Figure  5 to  be  typical.  Thus,  in  any  batch  of  thermistors  we  would  expect 
65-75%  of  them  to  have  drift  rates  less  than  5 m°  C/year  and  about  10%  of 
them  to  have  drift  rates  greater  than  10  m°  C/year.  A few  of  these 
would  have  drift  rates  high  enough  to  make  them  unuseable  for  any 
scientific  purpose.  We  should  point  out  that  high  drift  rates  seem  to 
decrease  with  time  but  we  will  need  several  more  years  of  calibrations 
before  we  can  elaborate  on  that  topic. 


DRIFT  RATE  MAGNITUDE  (m°C/YEAR) 

Figure  5 
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6.  Temperature  Accuracy  and  Precision 
6.  A.  NBS  Intercomparison 

In  December  1973  we  sent  17  thermistors  to  the  National  Bureau  of 
Standards  (NBS)  where  they  were  calibrated  by  Ms.  Sharrill  D.  Wood. 

Upon  their  return  to  W.H.O.I.  in  April  1974  they  were  calibrated  by 
Karl  Schleicher.  The  results  give  us  a comparison  of  the  results  of 
our  calibration  techniques  vs.  those  at  NBS. 

The  thermistors  chosen  for  this  intercomparison  were  of  two  types, 
both  preselected  by  the  manufacturer  for  stability  and  both  purchased 
for  a special  application  of  the  VACM.  One  type  ("5000  series")  has 
4000  ft  resistance  at  25°  C and  is  used  in  a differential  circuit  to 
measure  the  difference  in  water  temperature  between  the  top  and  bottom 
of  the  VACM,  thus  giving  local  temperature  gradient  (See  Section  6.  B., 

IWEX) . The  other  type  ("6000  series")  is  different  from  our  standard 
YSI  thermistor  in  that  it  is  pretested  for  stability  and  has  a slightly 
different  mechanical  configuration.  Their  temperature-resistance 
specifications  are  the  same  as  the  standard  VACM  thermistors,  having 
30,000  ft  resistance  at  25°  C. 

Normal  W.H.O.I.  calibrations  are  run  with  the  same  current  passed 
through  the  thermistor  as  in  the  VACM  at  that  temperature  so  that  the 
self-heating  will  be  quite  simixar.  For  this  inter comparison , however, 
all  resistance  measurements  at  NBS  and  W.H.O.I.  were  made  with  10  pa 
current  through  the  thermistors,  substantially  lower  than  normal  operating 
currents,  to  give  a minimum  amount  of  self-heating. 

Table  4 shows  the  difference  in  resistance  at  several  temperatures 
between  the  NBS  and  W.H.O.I.  calibrations  and  the  equivalent  temperature 
difference.  The  resistances  were  calculated  from  the  A,  B,  C constants 
(See  ApDendix  II)  at  the  indicated  temperatures.  The  temperature  difference 
was  calculated  by  multiplying  the  resistance  difference  by  dT/dR  as 
calculated  from  the  A,  B,  C constants  with  the  expression 

[B  + 3CUnR)2] 

dR  R 

The  individual  temperature  and  resistance  differences  in  Table  4 
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appear  fairly  random.  For  the  most  part  the  r.m.s.  differences  for 
each  thermistor  are  within  5 m°  C.  The  mean  difference  is  -1.6  ± 3.1  m°  C 
This  is  probably  what  one  should  realistically  expect  for  comparisons 
between  two  laboratories  and  may  be  indicative  of  the  absolute  accuracy 
of  calibrations  at  any  one  laboratory.  The  standard  deviation  is 
consistent  with  the  scatter  we  find  between  several  calibrations  of  a 
stable  thermistor. 

Calibration  results  of  two  of  the  "6000  series"  thermistors  were 
omitted  since  later  calibrations  showed  them  to  be  drifting  at  a rate 
which  would  make  comparison  of  the  NBS  and  W.H.O.I.  results  meaningless. 
Four  of  the  thermistors  listed  were  destroyed  for  testing  purposes 
after  only  two  calibrations.  These  were  6003,  6017,  6046,  6118. 

These  show  some  of  the  highest  differences  for  the  6000  series  thermistors 
in  Table  4.  Because  of  the  limited  number  of  calibrations  we  cannot 
be  sure  that  the  differences  were  not  caused  by  drift. 


. - 


Aik.  * 


26 


6 . B . IWEX 

During  November  1973  a stable  tripod  mooring  was  deployed  by  the 
Buoy  Project  at  27°  43'N,  69°  51'W  in  5455  m water  depth.  The  mooring  was 
in  place  for  44  days.  On  the  IWEX  (for  Internal  Wave  Experiment)  mooring 
were  mounted  17  VACMs,  9 temperature/pressure  recorders,  and  4 model  850 
current  meters.  The  VACMs  were  modified  to  measure  temperature  difference 
between  thermistors  mounted  in  pods  at  the  top  and  bottom  of  each  instru- 
ment as  well  as  absolute  temperature.  Because  of  the  temperature  difference 
measurements  extensive  calibration  operations  were  performed.  These  calib- 
ration operations,  combined  with  the  stability  of  the  mooring  and  the  close 
spacing  of  a cluster  of  7 VACMs  at  the  top,  offered  an  excellent  chance  to 
determine  the  achievable  accuracy  in  the  VACM  measurement  of  absolute 
temperature  in  situ. 

Figure  6 illustrates  the  configuration  of  the  IWEX  mooring.  The 
depths  of  the  levels  and  the  horizontal  separations  of  instruments  on  two 
different  legs  but  at  the  same  level  are  given  in  Table  5 as  well  as  the 
temperature  gradient  at  each  level  derived  from  a mean  temperature  profile 
(Millard,  1974;  Tarbell,  Briscoe  and  Chausse,  1975)  . 

Calibrations 

The  thermistors  used  in  the  absolute  temperature  circuit  of  the 
IWEX  VACMs  were  manufactured  by  Thermometries,  Inc.  Mechanically  and 
electrically  they  are  quite  similar  to  the  YSI  thermistors  we  had  used 
previously.  Instead  of  being  mounted  within  the  current  meter,  they  were 
mounted  externally  in  pods  to  decrease  the  time  constant. 

Before  delivery  the  thermistors  were  stabilized  by  temperature  cycling 
and  then  tested  for  maximum  stability  by  repeated  calibrations  at  regular 
intervals  by  the  manufacturer.  In  addition,  they  were  calibrated  in  our 
own  facility  (see  Section  4 of  this  report)  before  and  after  the  experi- 
ment and  twice  more  in  the  two  years  since.  We  have  not  included  these  in 
the  thermistors  reported  on  in  Section  5 because  of  this  limited  calibration 
history. 

Because  of  the  stringent  requirements  of  the  temperature  difference 
measurements  the  thermistor  pods  were  immersed,  several  at  a time,  in  a 
constant  temperature  bath  before  and  after  the  experiment  for  system  calib- 
ration. The  bath  was  run  at  0,  15,  30°  C while  the  thermistors  remained 


connected  to  the  current  meters.  These  data  were  recorded  on  the  IWEX 
data  tape  for  each  instrument. 

The  constant  temperature  bath  was  a Tronac  Model  400,  modified  by 
the  addition  of  one  inch  of  urethane  insulation  on  all  sides.  The  tem- 
perature controller  was  a Tronac  Model  40  used  in  conjunction  with  an 
external  cooling  unit.  These  provided  excellent  temperature  control  for 
the  calibrations,  maintaining  control  to  within  2 m°  C (averaged  over 
225  seconds)  for  periods  of  hours  with  temperature  differences  through 
the  tank  of  less  than  .2  m°  C.  Absolute  temperature  of  the  tank  was 
measured  with  a platinum  resistance  thermometer  with  calibrations  trace- 
able to  the  National  Bureau  of  Standards. 

The  noise  level  and  short-term  stability  of  the  instruments,  and 
the  stability  of  the  calibration  bath  were  within  design  specifications. 

As  a test  five  instruments  were  run  overnight  in  the  bath.  The  absolute 
temperatures  recorded  by  the  current  meters  were  constant  within  ±1  data 
count  (±.2  m°  C)  averaged  over  225  seconds. 

Results 

Biases  and  drift  rates  of  the  thermistors  were  computed  from  the 
bath  tests  and  our  standard  thermistor  calibrations.  Bias  is  the  error 
in  measurement  of  absolute  temperature  by  a current  meter  and  thermistor 
at  the  time  of  the  pre-mooring  bath  test.  Drift  rate  is  computed  from  the 
difference  in  error  of  measurement  between  the  pre-mooring  and  post-mooring 
bath  tests  and  the  pre-mooring  and  post-mooring  thermistor  calibrations. 
Only  two  of  the  seventeen  thermistors,  those  at  the  positions  B2  and  C5, 
showed  significant  drift  and  needed  to  be  corrected  but  the  magnitudes 
of  these  two  drift  rates,  9 and  11  m°  C/month,  were  surprisingly  large 
in  view  of  the  pretesting  performed  by  the  manufacturer. 

Table  6 lists  means  of  absolute  temperature  for  all  the  levels  where 
there  was  more  than  one  current  meter.  They  are  listed  as  ten  consecutive 
4-day  means  and  the  overall  40-day  mean.  The  B2  and  C5  temperatures  have 
been  corrected  for  bias  and  drift.  No  corrections  were  required  for  the 
rest  of  the  temperatures. 
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In  Table  7 are  the  temperature  differences  at  each  level  for  the 
4-day  means,  the  40-day  meah,  and  the  40-day  mean  divided  by  the  tem- 
perature gradient  at  that  level  from  Table  5.  If  we  assumed  that  the 
temperature  sensors  were  all  perfect  then  this  last  quantity  would  be  the 
mean  depth  difference  between  sensors  on  one  level. 

Now  look  at  Table  7 in  detail  keeping  in  mind  our  stated  temperature 
accuracy  of  10  m°  C.  Sensors  A1  and  Cl  are  only  6.1  meters  apart  and  the 
temperatures  agree  within,  at  most,  2 m°  C.  At  level  2,  with  a separation 
of  8.5  meters,  A2  and  C2  agree  to  no  worse  than  3 m°  C.  B2  is  one  of  the 
two  thermistors  which  had  a rather  large  drift  rate  and  is,  therefore, 
not  as  reliable  as  the  others.  In  spite  of  this,  it  agrees  with  A2  and  C2 
within  5 m°  C except  for  one  of  the  4-day  periods.  A4  and  B4  agree  within 
5 m°  C.  At  levels  5,  6,  10  with  substantially  larger  sensor  separations 
the  B leg  temperatures  are  substantially  higher  than  for  legs  A or  C. 

To  explain  this  requires  that  the  sensors  on  the  B leg  be  1,  2,  3 meters 
higher  for  the  5,  6,  10  levels  than  the  sensors  on  the  other  two  legs. 
Although  there  were  not  enough  working  pressure  sensors  on  the  mooring  to 
define  the  motions  we  can  rationalize  the  direction  if  not  the  magnitude, 
of  the  leg  motions. 

Anchors  at  the  base  of  legs  A and  C were  lined  up  along  a north- 
south  line.  The  currents  in  the  top  1500  m were  predominantly  from  west 
to  east  during  the  whole  40  days.  This  would  cause  legs  A and  C to  be 
depressed  relative  to  leg  B.  The  current  veered  somewhat  about  the  easterly 
direction  with  time,  however,  so  we  cannot  say  much  about  the  relative 
depths  of  sensors  on  legs  A and  C.  It  is  apparent  from  Table  7 that  the 
B-A  and  B-C  temperature  differences  are  substantially  larger  than  the  C-A 
at  levels  5,  6,  and  10. 

1/  The  absolute  temperature  data  from  the  IWEX  mooring  are  consistent 

with  our  claimed  temperature  accuracy  of  10  m°  C.  They  are  particularly 
supportive  where  the  horizontal  separations  are  small  and  we  can  expect 
the  vertical  separations  to  be  less  than  about  50  cm. 


IWEX  Mooring  Data 


Level 

Depth 

(m) 

Horizontal 

Separation 

of 

Instruments 

(m) 

Temperature 
Gradient 
(m°  C/m) 

1 

603.6 

6.1 

18.6 

2 

605.7 

8.5 

18.6 

4 

610.6 

14.0 

18.8 

5 

639.5 

44.0 

19.8 

6 

730.6 

139 

22.5 

8 

1014.4 

441 

12.2 

10 

1023.1 

450 

11.7 

14 

2050.4 

1600 

1.4 

Table  6 
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IWEX  4-Day  Mean  Temperatures,  B2  and  C 5 Corrected  for  Thermistor  Drift 


40-Day 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Mean 

A1 

14.029 

13.964 

13.968 

14.015 

14.116 

14.158 

14.180 

14.186 

14.181 

14.173 

14.097 

Cl 

14.027 

13.962 

13.967 

14.014 

14.116 

14.157 

14.179 

14.184 

14.180 

14.172 

14.096 

A2 

13.995 

13.927 

13.930 

13.975 

14.079 

14.120 

14.144 

14.152 

14.144 

14.141 

14.061 

B2 

13.988 

13.922 

13.925 

13.972 

14.075 

14.117 

14.141 

14.150 

14.141 

14.139 

14.057 

C2 

13.995 

13.926 

13.930 

13.974 

14.078 

14.119 

14.142 

14.149 

14.141 

14.138 

14.059 

A4 

13.865 

13.798 

13.793 

13.836 

13.946 

13.985 

14.008 

14.024 

14.005 

14.015 

13.927 

B4 

13.866 

13.800 

13.794 

13.839 

13.949 

13.990 

14.011 

14.027 

14.009 

14.017 

13.930 

A5 

13.288 

13.231 

13.181 

13.241 

13.358 

13.429 

13.452 

13.461 

13.429 

13.463 

13.353 

B5 

13. 301 

13.249 

13.198 

13.262 

13.380 

13.453 

13.473 

13.484 

13.448 

13.482 

13.373 

C5 

13.289 

13.237 

13.188 

13.250 

13.367 

13.436 

13.457 

13.467 

13.430 

13.468 

13. 359 

A6 

11.254 

11.155 

11.142 

11.134 

11.264 

11.389 

11.383 

11.413 

11.431 

11.442 

11.301 

B6 

11.293 

11.198 

11.184 

11.183 

11.317 

11.447 

11.435 

11.464 

11.488 

11.496 

11.350 

C6 

11.260 

11.162 

11.150 

11.147 

11.276 

11.387 

11.371 

11.402 

11.410 

11.434 

11.300 

A10 

6.196 

6.219 

6.244 

6.269 

6.271 

6.317 

6.437 

6.594 

6.617 

6.583 

6.375 

BIO 

6.219 

6.247 

6.278 

6.303 

6.315 

6.365 

6.476 

6.635 

6.669 

6.627 

6.413 

CIO 

6.192 

6.222 

6.252 

6.280 

6.278 

6.312 

6.422 

6.575 

6.594 

6.572 

6.370 

A14 

3.538 

3.563 

3.605 

3.605 

3.580 

3.581 

3.620 

3.635 

3.629 

3.619 

3.597 

B14 

3.554 

3.589 

3.628 

3.629 

3.615 

3.602 

3.641 

3.664 

3.655 

3.639 

3.622 

C14 

3.547 

3.597 

3.622 

3.621 

3.598 

3.592 

3.626 

3.641 

3.626 

3.619 

3.609 

Table  7 


1 

2 

IWEX  Temperature  Differences 
3 4 5 6 7 8 9 

10 

Means 

Equivalent 

Depth 

Cl-Al 

2 

2 

1 

1 

0 

1 

1 

2 

1 

1 

1 

Differences 

(m) 

.1 

C2-A2 

0 

-1 

0 

-1 

-1 

-1 

-2 

-3 

-3 

-3 

-2 

.1 

B2-A2 

-7 

-5 

-5 

-3 

-4 

-3 

-3 

-2 

-3 

-4 

-4 

.2 

B2-C2 

-7 

-4 

-5 

-2 

-3 

-2 

-1 

-1 

0 

-1 

-2 

.1 

B4-A4 

1 

2 

1 

3 

3 

5 

3 

3 

4 

2 

3 

-.2 

C5-A5 

1 

6 

7 

9 

9 

7 

5 

6 

1 

5 

6 

-.3 

B5-A5 

13 

18 

17 

21 

22 

24 

21 

23 

19 

19 

20 

-1.0 

B5-C5 

12 

12 

10 

12 

13 

17 

16 

17 

18 

14 

14 

-.7 

C6-A6 

6 

7 

8 

13 

12 

-2 

-12 

-11 

-21 

-8 

-1 

0 

B6-A6 

29 

43 

42 

49 

53 

58 

52 

51 

57 

54 

49 

-2.2 

B6-C6 

33 

36 

34 

36 

41 

60 

64 

62 

78 

62 

50 

-2.2 

C10-A10 

-4 

+3 

+8 

11 

7 

-5 

-15 

-19 

-23 

-11 

-5 

.4 

B10-A10 

23 

28 

34 

34 

44 

48 

39 

41 

52 

44 

38 

-3.2 

BIO-CIO 

27 

25 

26 

23 

37 

53 

54 

60 

75 

55 

43 

-3.7 

C14-A14 

9 

34 

17 

16 

18 

11 

6 

6 

-3 

0 

12 

8.6 

B14-A14 

16 

26 

23 

24 

35 

21 

21 

29 

26 

20 

25 

17.9 

B14-C14 

7 

-8 

6 

8 

17 

10 

15 

23 

29 

20 

13 

9.3 

I 
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TRIMOOR  GEOMETRY 


6.  C.  Mooring  551  Intercomparison 

In  December  1974,  a mooring  was  implanted  at  Site  J,  36°N,  70°W, 
with  four  VACMs  rigidly  shackled  together  at  a nominal  2000  m depth. 
Mooring  551  was  in  place  for  128  days  during  which  all  four  instruments 
worked  satisfactorily.  A fifth  instrument,  containing  a pressure  sensor, 
did  not  work.  Because  of  this,  the  depths  of  the  VACMs  are  known  only 
to  about  ±20  m. 

Since  the  mooring's  purpose  was  a detailed  comparison  of  the  VACMs, 
the  depth  of  the  instrument  grouping  was  chosen  to  be  in  a region  of  low 
gradients.  Figure  7 shows  temperature  vs.  depth  from  a CTD  station  made 
at  the  mooring  location  when  it  was  set.  The  inset  shows  a more  detailed 
plot  for  the  depth  range  1900-2100  m,  a straight  line  least  squares  fitted 
to  the  CTD  data,  and  four  points  representing  the  128  day  temperature 
means  for  the  four  VACMs.  The  straight  line  has  a slope  of  -0.71  m°  C/m, 
or  the  equivalent  of  a 4.3  m°  C spread  in  temperature  over  the  6 m between 
the  top  and  bottom  thermistors. 

Table  8 shows  the  128  day  temperature  means  from  the  four  VACMs 
and  their  standard  deviations.  For  each  7.5  minute  recording  interval  a 
mean  temperature  was  computed  from  the  four  individual  temperatures  and 
the  difference  between  the  mean  and  the  individual  temperatures  computed 
to  see  how  well  the  four  instruments  tracked.  Row  4 is  the  128  day  time 
mean  of  these  differences  for  each  instrument  and  row  5 is  the  standard 
deviation  of  this  quantity.  The  results  show  that  the  instruments  track 
very  well  and  that  the  agreement  of  the  means  is  not  fortuitous. 


Table  8 

Mooring  551  Temperatures  Averaged  over  128  Days 


1. 

Instrument 

5512 

5513 

5514 

5515 

2. 

Mean  temperature 

3.935°  C 

3.933°  C 

3.936°  C 

3.932° 

3. 

Standard  deviation 

0.024 

0.024 

0.024 

0.024 

4. 

Mean  difference 

0.001 

-0.001 

0.002 

-0.002 

5. 

Standard  deviation 

0.002 

0.001 

0.001 

0.002 

It  would  be  appealing  to  attempt  some  kind  of  a direct  comparison 
between  absolute  temperatures  from  VACMs  and  CTD  but  the  uncertainty 
in  current  meter  depth  rules  this  out.  The  CTD  data  allow  us  to  say 
that  the  time  averaged  temperature  gradient  in  the  vicinity  of  2000  m is 
probably  very  close  to  -0.7  m°  C/m.  The  absolute  temperatures  are 
consistent  with  the  CTD  temperatures  to  within  the  accuracy  of  the  VACM 
depths . 
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6.  D.  Error  Analysis 

In  the  previous  sections  of  this  report  we  have  described  the 
components  of  the  VACM  temperature  measuring  system,  some  of  the  errors 
in  those  systems,  and  intercomparisons  which  have  indicated  some  of  the 
relative  uncertainties  in  our  temperature  measurements.  We  have  not 
been  able  to  make  an  in  situ  comparison  with  a different,  and  preferably 
certifiably  accurate,  instrument,  nor  have  we  been  able  to  test  the 
absolute  accuracy  in  the  laboratory.  In  this  section  we  will  compute 
analytically  the  accuracies  of  both  relative  and  absolute  temperature 
measurements.  To  the  extent  that  our  computed  accuracy  of  relative 
measurements  agrees  with  the  intercomparison  results  we  may  have  some 
confidence  in  the  computed  absolute  accuracy. 

In  Table  9 is  a list  of  possible  sources  of  error  with  an  estimate 
of  the  contribution  of  each  to  the  instrument  bias  in  the  measurements 
and  the  measurement  variation. 

Tabla  9 

Error  Sources  in  VACM  Temperature  Measurements 

Instrument 
Bias 

Thermistor 

a.  Calibration  and  drift  rate 

b.  Linearity  +.3m°C 

c.  Self  heating  +lm°C 

V/F  Converter 

a.  Stability 

b.  Linearity 
Rest  of  VACM 

a.  Time  base  stability  (900  second  record- 

ing interval) 

b.  Integrator 


Measurement 

Variation 


±4m°C 

0 

0 

lm°C 

lm°C 

±.3m°C 

±.06m°C 


4|*. 
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The  following  describes  the  terms  in  the  table  in  more  detail  and 
justifies  the  uncertainties  given. 

Thermistor 

a.  Calibration  and  drift  rate  - From  the  NBS  intercomparison  we 
feel  sure  that  the  absolute  accuracy  of  each  of  our  calibrations  is 
about  ±3  m°  C.  The  scatter  about  the  line  fitting  calibration  vs. 
time  data  is  ±1  m°  C or  less.  The  total  uncertainty  in  expressing  T 
as  a function  of  R by  the  equation 

T_ 1 = A + B(dnR)  + C (inR) 3 


i 


<1 

< 


is  then  about  ±4  m°  C. 

b.  Linearity  - Thermistors  are  nonlinear  sensors.  Computing  a mean 
temperature  from  a mean  thermistor  resistance  introduces  an  error  which 
depends  on  the  amplitude  of  fluctuations  during  the  averaging  interval. 
Typical  temperature  records,  both  in  and  below  the  thermocline  show 
temperature  standard  deviations  of  about  200  m°  C.  For  a sine  wave 
fluctuation  of  this  amplitude  the  error  would  be  + .3  m°  C.  This  is  a 
systematic  and  not  a random  error  so  it  contributes  to  the  instrument  bias. 

c.  Self  heating  - In  the  input  circuit  of  the  V/F  converter,  a 
precise  3.9  volts  is  impressed  across  the  thermistor  and  a precision 
37.3  KQ  fixed  resistor  is  connected  in  series  with  the  thermistor.  In 
Table  10  are  shown: 

- nominal  temperature 

- mean  resistance  and  standard  deviation  for  the  thermistors 
in  Table  3 with  drift  rates  less  than  3.5  m°  C/year. 

- power  dissipated  in  thermistor  with  standard  deviation 
due  to  deviation  in  R. 

- thermistor  self-heating  with  assumed  dissipation  constant  of 
3 mw/°C  and  standard  deviation  due  to  deviation  in  R. 

- volts  across  thermistor  resistance  R in  current  meter 
iiircuit  with  standard  deviation  due  to  deviation  in  R. 

- volts  applied  across  thermistor  in  calibration  equipment. 


T 

R 


th 


V 


calib 


dT/dV  - sensitivity  of  indicated  temperature  to  an  error  in  the 


voltage  applied  across  the  thermistor. 


Table  10 

Thermistor  Self-Heating 
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T 

R(ohms) 

p (yw) 

At  (m°  C) 

V . (volts 
th 

V (volts 

cal 

dT 

dV 

m°  C, 
volt' 

0°C 

95080±377 

82.521.14 

27.511.05 

2.8011.003 

2.80 

20 

5 

74453±286 

90.681.12 

30.231.04 

2.5981.003 

2.60 

23 

10 

58712+220 

96.871.08 

32.291.03 

2.3851.003 

2.39 

27 

15 

46611+171 

100.691.04 

33.561.01 

2.166+.004 

2.17 

31 

20 

37244+135 

101.941.0003 

33.981.0001 

1.9491.004 

1.95 

35 

25 

299471106 

100.731.04 

33.581.01 

1.7371.003 

1.74 

39 

30 

24222184 

97.341.07 

32.451.02 

1.5351.003 

1.54 

42 

Since  the  same  voltage  is  applied  to  the  thermistor  during  calibration  as 
in  the  current  meter,  the  bias  due  to  self-heating  is  very  nearly  the 
same  in  the  two  situations.  This  net  bias,  therefore,  does  not  appear 
in  Table  9.  There  are  three  kinds  of  errors  that  can  occur  through  the 
self-heating: 

1.  An  error  in  the  volts  applied  to  the  thermistor  during  calibration 
can  cause  an  amount  of  self-heating  differing  from  that  in  the  current 
meter.  An  error  of  0.1  volt  would  cause  up  to  5 m°  C error,  indicating 
that  a modicum  of  care  is  required. 

2.  Because  various  thermistors  have  somewhat  different  resistances 

at  the  same  temperature,  the  voltage  across  them,  and  thus  the  self-heating 
varies.  For  the  range  of  thermistors  in  this  study,  however,  the  difference 
is  at  most  .05  m°  C,  a negligible  amount. 

3.  The  assumed  dissipation  constant  of  3 mw/°C  is  a nominal  value. 

The  important  thing  is  to  have  the  dissipation  constant  the  same  in 

the  current  meter  and  in  the  calibration  setup.  We  have  not  been  able 
to  measure  either  dissipation  constant  but  have  tried  to  make  the  sit- 
uations quite  similar  physically.  We  feel  confident  that  the  error  due 
to  difference  in  self-heating  is  less  than  ±1  m°  C. 

V/F  Converter 

a.  Stability  - We  saw  in  Section  3 that  the  V/F  converter  is  stable 
to  3 parts  in  105,  equivalent  to  about  ±1  m°  C. 

b.  Linearity  - We  saw,  also  in  Section  3,  that  the  period  of  the 
V/F  converter  output  signal  is  a linear  function  of  thermistor  resistance 
to  within  the  equivalent  of  ±1  m°  C. 


~ iWi^a«i  Filin  i Vi  ill 
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Rest  of  VACM 

a.  Time  base  stability.  The  stability  of  the  crystal  oscillator 
which  provides  the  time  base  in  the  VACM  is  stable  to  1 part  in  10s, 
equivalent  to  .009  sec  for  a recording  interval  of  900  sec.  The  un- 
certainty in  temperature  due  to  this  uncertainty  can  be  derived  from 
the  three  equations: 


T = [A  + B(fcnR)  + C(fcnR) *1_1 


<?  ♦ V Rs 

K:  - k2  - i 


p - i 

N 

where : 

A,  B,  C - thermistor  calibration  constants 
R - thermistor  resistance  at  temperature  T 
= fixed  resistor,  see  section  2 

= V/F  converter  constants,  see  section  2 
t = current  meter  recording  interval 

N = counts  accumulated  during  T by  counter  on  output  of  the  V/F 
converter 

The  uncertainty  in  temperature  due  to  variations  in  the  time  base  is  given  by: 

at|  At 

!n  9r  9p  9t 


" T [B  + 3c(£nR)  1 

9r  Rs(K2  - V 

9p  [p(k  - k2)  - l]2 


9p  = l 

9t  N 
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Typical  values  of  the  variables  at  10°  C are: 

24270  £2 
94980  fi 

1.0036  x lo'J  sec 
2.1548  x 10-5  sec 
.009  sec 

Substituting  these 

At | = ±.3  m°  C 

N 


T = 283.15°  C 
R = 58712  fi 
B = .22519048  x lo" 
C = .114084  x 10-6 
N = 575117 


R.  = 


R„  = 


P_  = 


2 

At 


b.  Integrator  - The  integrator  is  just  a counter  and  has  an 

uncertainty  of  ±1  count.  From  Table  1 we  see  that  this  is  equivalent 

to  ±.06  m°  C,  a negligibly  small  uncertainty. 

The  total  variance  of  temperature  is  the  sum  of  the  squares  of  the 

uncertainties.  For  the  measurement  variations, 

Var  = 18 (m°  C) 2 
T 

The  relative  uncertainty  we  expect  in  the  temperature  measurements  is 
then  of  the  order  of  ±4.2  m°  C which  agrees  well  with  the  results  of 
the  IWEX  inter comparison. 

Adding  the  contributions  of  instrument  biases  to  the  variance, 


VarT  = 19 (m°  C) 2 

the  uncertainty  we  expect  in  the  absolute  temperature  measurements  is 
then  about  ±4.4  m°  C. 

We  feel  that  this  uncertainty  estimate  is  a little  optimistic  but 
not  completely  out  of  line.  Nelson  Hogg  (personal  communication)  after 
contouring,  separately,  temperatures  recorded  during  MODE  by  current 
meters  and  CTDs,  stated  that  the  two  sets  appeared  to  agree  within  ±10  M° 
Our  best  guess  as  to  the  achievable  accuracy  of  absolute  temperatures 
with  the  VACM  then,  is  perhaps  better  than  ±10  m°  C but  no  better  than 
±4.5  m°  C. 


■ ^ * i 


41 


References 


Bennett,  A.  S.,  1972.  The  calibration  of  thermistors  over  the  temperature 
range  0-30°  C.  Deep-Sea  Res. , 19,  157-63. 

Bradshaw,  A. , and  K.  E.  Schleicher,  1970.  Direct  measurement  of  thermal 
expansion  of  sea  water  under  pressure.  Deep-Sea  Res . , 17,  691-706. 

Millard,  R. , Jr.,  1974.  CTD  averages  from  IWEX.  MODE  Hot  Line  News, 
no.  52,  3.  (unpublished  manuscript). 

Stei.ihart,  J.  S.,  and  S.  R.  Hart,  1968.  Calibration  curves  for  thermistors. 
Deep-Sea  Res.,  15,  497-503. 

Tarbell,  S.,  M.  G.  Briscoe,  and  D.  Chausse,  1975.  A compilation  of  moored 
current  data  and  associated  oceanographic  observations.  Volume  IX 
(1973  Internal  Wave  Experiment  (IWEX)).  W.H.O.I.  Ref.  75-68 
(unpublished  manuscript) . 

Vector  Averaging  Current  Meter  Technical  Manual,  AMF  Publication  #SLS  106-11419. 
AMF  Electrical  Products,  Development  Division,  1025  North  Royal  Street, 
Alexandria  VA. 


42 


Appendix  I 
Some  Useful  Numbers 

Some  useful  information  is  summarized  in  Table  1-1. 

Column  1 - The  seven  nominal  temperatures  at  which  we  monitor  ther- 
mistor and  V/F  converter  characteristics. 

Column  2 - Nominal  resistances  at  these  temperatures  of  the  YSI 
44032  thermistor. 

Column  3 - Nominal  dT/dR  for  YSI  44032  thermistors. 

Column  4 - Nominal  period  of  the  V/F  converter  output  signal  for 
the  Column  2 resistances.  Our  V/Fs  are  normally  adjusted  to  agree  with 
the  nominal  periods  within  ±.01  ys  at  74.440Q. 

Column  5 - Some  typical  VACM  recording  intervals. 

Column  6 - Numbers  of  counts  which  would  be  accumulated  during  the 
various  recording  intervals  with  thermistor  resistances  of  Column  2. 

Column  7 - Temperature  resolution  due  to  1 count  resolution  in  the 
total  counts. 

Column  8 - Uncertainty  in  temperature  due  to  uncertainty  in  either 
PI  or  P2  of  +1  ys. 

Column  9 - Accuracy  requirement  for  PI,  P2  for  a ±1  m°  C accuracy 
in  temperature. 

Columns  10,  11,  12  - Accuracy  requirements  in  A,  B,  C thermistor 
constants  for  ±1  m°  C temperature  accuracy.  See  Section  5 for  definition 
of  A,  B,  C. 
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Appendix  II 
CALIBB  Program 

Included  in  this  appendix  are  a description  of  the  CALIBB  program, 
a listing  of  the  program,  and  an  example  of  input  and  output  data  from 
a calibration.  CALIBB  was  written  using  Xerox  FORTRAN  IV  and  contains 
some  statements  which  are  not  compatible  with  other  versions  of  FORTRAN  IV. 
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CALIBB, 1 


NAME:  CALIBB 

TYPE:  Main  program  - processor 

PURPOSE : To  process  thermistor  calibration  data 

MACHINE:  XDS  Sigma  7 

SOURCE  LANGUAGE:  Xerox  FORTRAN  IV 

PROGRAM  CATEGORY:  Utility 

DESCRIPTION: 

CALIBB  was  written  to  fit  thermistor  calibration  data  to  the 
equation 


1/T  = A + B *LN (R)  + C* (LN (R) ) **3 

where  T is  absolute  temperature  in  °K  and  R is  thermistor 
resistance . 

CALIBB  provides  a choice  of  three  inputs; 

1.  Direct  values  of  temperature  and  corresponding  thermistor 
resistance. 

2.  Output  of  a Guildline  bridge  connected  to  a platinum  resistance 
thermometer  is  converted  to  temperature.  Thermistor  resistance 
has  corrections  applied  for  lead  resistance  and  bridge  dial 
resistance. 

3.  Output  of  a Mueller  bridge  connected  to  a platinum  resistance 
thermometer  is  converted  to  temperature. 

The  result  of  all  three  types  of  input  data  is  a set  of  temperatures 
and  thermistor  resistances  for  each  thermistor.  These  are  then 
fitted  to  the  above  equation  by  the  method  of  least  squares.  The 
resulting  A,  B,  C constants  are  used  for  computing  temperature 
from  thermistor  resistance. 

CALIBB  was  written  particularly  to  process  the  calibration  data 
produced  by  K.  Schleicher  and  A.  Bradshaw  for  Buoy  Project  ther- 
mistors used  in  current  meters. 

INPUT:  By  cards 

For  each  batch 

Card  1 is  the  same  for  all  three  types  of  input  data. 

1 card  - Batch  number,  type  of  input  (DIRECT,  GUILD, MUELL) , number 
of  thermistors  in  batch,  number  of  nominal  temperatures 
per  thermistor,  date  of  calibration,  source  of  calibration  data 

(T7,I3,  T16,A4,  T36.12,  T51,12,  2X,A2,  IX, A3,  IX, A2,  1X,3A4) 


..4 
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CALIBB , 2 


Type  1,  Direct 

For  each  thermistor 

1 card  - Header  card  - thermistor  number,  batch  number 
(T13 ,14 , T25 , 13) 

1 or  more  cards  - Data  - 4 pairs  of  temperature , resistance 
(4  (2F) ) 

Type  2,  Guildline  bridge 
For  each  batch 

1 or  more  cards  - list  of  thermistor  number  in  calibration 
batch  in  order  of  subsequent  data  appearance 
(161) 

1 card  - Platinum  probe  identifier,  zero  ratio  of  bridge 
(T9,A4,  T25 , F) 

For  each  nominal  temperature 

1 card  - nominal  temperature,  data  date,  initial  bridge 
ratio,  volts  across  thermistor,  number  of  book  containing 
original  data,  page  in  book 

(T9,F4. 1 , T14 , 3A4 , T36,F,  T52,F,  T63I,  T74,I) 

1 or  more  cards  - data  pairs  consisting  of  position  in  bath, 
bridge  ratio,  and  uncorrected  resistance 
(T8 , 12 , T10 , 2F ) 

Type  3,  Mueller  bridge 

For  each  batch 

1 or  more  cards  - list  of  thermistor  numbers  in  calibration 
batch  in  order  of  subsequent  data  appearance 

(161) 

For  each  nominal  temperature 

2 cards  - nominal  temperature,  data  date,  zero  error  of  bridge, 
volts  across  thermistor,  book  containing  original  data,  page 

in  book 

(T9  ,F  , T13 , 3A4 , T34,F,  T52,F,  T63,I,  T74,I) 

1 card  - NRML1MA, RVRS1MA ,RVRS14MA ,COEFF , RZERO  (For  definitions 
see  program  listing) 

(T9  ,F  , T26 ,F  , T44,F,  T59,F,  T73,F' 

For  each  thermistor 

1 card  - data  for  this  nominal  temperature;  BRIJN1MA,NULLN1MA, 
BRIJR1MA , NULLR1MA , BRIJM14MA ,NULLM14MA , RES 

(T9 , 7F ) 
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CALIBB, 3 


OUTPUT: 

For  GUILDLINE  or  MUELLER  input 

Printer:  Listing  of  computed  platinum  resistance  temperatures  and 

corrected  thermistor  resistances  by  nominal  temperature 
for  data  validation 

Card  punch:  Temperatures  and  resistances  for  each  thermistor 
for  data  archiving 

For  all  three  input  types 

Printer:  Listing  of  A,B,C  constants  and  statistics  of  fit  for 

each  thermistor 

Card  punch:  Three  duplicate  cards  for  each  thermistor  containing 

thermistor  number,  batch  number,  calibration  data,  A,B,C 
constants,  temperature  RMS  residual  (measure  of  fit) 


USEAGE : 

For  processing  using  the  source  deck 

! LIMIT  (TIME, 2) , (CORE, 20) 

(MESSAGE  PUNCHES  CARDS 
! INTERP  FULL 
(FORTRAN  LS , GO 
(LOAD  (GO),  (UNSAT, (3)) 

SOURCE  DECK 

(RUN 

(DATA 

DATA  CARDS 
RESTRICTIONS:  None 

STORAGE  REQUIREMENTS:  Not  applicable 

SUBPROGRAMS  REQUIRED:  None 

OPERATIONAL  ENVIRONMENT:  CP-V  Monitor 

TIMING: 


Batch  of  24  thermistors  and  7 nominal  temperatures  requires  about 
3.9  charge  units  when  run  with  the  source  deck,  2.2  charge  units 
when  run  with  an  object  deck 
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PR03RAH  CALtBB 
BY  R.  payne  £T  al 
8 SEPTEMBER  t976 


Calis  pits  a set  of  twehmistor  calibration  data  to  t«e  equation 
1/t.a*h*UN(W)*C*(LN(«>**3> . CALIB  PERFORMS  Twr  FOLLOW I NG  TASKS! 

1.  WE ADS  TWE  INPUT  DATA 

?.  computes  temperatures  from  t«e  platinum  resistance  thermometer  <prtj 
bridge  Headings  and  prints  but  tme  Resui  tins  r/t  values  for  Data 

validation.  PROVISION  is  MADr  FOR  INPuT  op  r,t  DATA  DIRECTLY  rOR  TASK  3 
3.  fits  T*E  P/T  Data  tB  TWE  ABOVE  Enu*TrON  AND  OUTPUTS  TWE  PESJ-TS  of 
the  Fit  F9R  EACH  INDIVIDUAL  THeRM^TOR, 

ordinarily.  Calibration  measurements  are  made  on  thermistors  at  seven 
TempEwATURPS,  nominally  0.30  C AT  5 c INTERVAL*. 

bOTP  a MufcLLpR  AND  a GUILULINe  bRIORe  HAve  “pen  USrD  TO  measURe  TW£  ®RT 

Resistance  although  only  the  guiloline  is  eD  at  present,  starting  with 
batch  5i;  corrections  are  applied  within  t«e  program  eor  lead  Resistance  and 

HRJDGe  DeCaDe  RESISTANCES. 


OUTLINE  t)F  CALIHB  and  SUBROUTINES 


*Rp  aD  Batch  number  (BatChnum),  input  SOURCe  type  (Qmdj,  NUMBeP  bp  tbp^mistORS 
jn  BatCh  inUMthfRMj,  number  »f  nomjnaL  temperaturES  ( nUMnsMteMP, 
calibration  date  (Calibdatej#  plaCe  bf  data  origin  (CalibplaCej . 

• TEST  FOR  E(*'3»FND  OF  DATA 

• TFSt  For  RFA5ftNAB|.F  VALUES  OF  NUMTwERM,  NUMNOMTEMP 

• TEST  GMO  EOR  GUILDLINE/  OR  mueLLER  BRIDGE,  OR  D I RrCT  R, T INPUT 

•Call  guild, muell,or  direct 
•EqR  F*Ch  TmFrM!sT0R 
•call  OBCFIT 
•Call  printout 

print  output  of  least  suuaRes  fit  for  bne  THpRMjST8R 
•call  punchout 

pUnCh  R*T  DATA  P8InTS  *nD  RESULtS  OF  lEAST  5OUARE5  EjT  E jr  EACh  TjErMisT8R 

• END 

subroutines  for  input  and  conversion  oe  bridge  data 
guildline  bridge 
•subroutine  guild 

•read  therhistor  bath  seuuencf  numbers  (thErmjd) 

•READ  °TPROBe  ■ WHICH  PLATINUM  PROBE  USED,  BR AD.BR AOSWAW,  BUOY.BJOY  LABj 

and  probe  bridge  ratio  at  „ c and  x m4  (ZeRoratioj 
•for  e aCw  nominal  tf.mperature 

•read*  nominal  temp  (n»mtemP), date, bridge  ratio  when  probe  is  at 
measurement  te^eRatureiinttRatiO). volts  across  thermist9R< volts>, 

BOOK  AND  PAGE  BE  DATA  JW  SCwLEICHeriS  NOTEBOOK, 

•EOR  e*CH  THERMISTOR 

*rE*D’  BATh  sEUuEnCE  numBEr(BAThIOuOIFFErEnCE  BETwECn  InITrATio 
and  actual  bridge  ratio(DeltRatio, uncorreCted  resistance  of 
therm ; s TO R(URES) • 


V 


#TEsT  Few  B ATCHNUM.  GT  « 51  r0 

• YES  3 

•CALL  RESCORR 

CORRECT  All  ThErM^Tor  REsISt*nCEs  F9r  LE*d  AmO  Bp  1 DUE  DIAL 
RESISTANCE 

• SET  PLATINUM  pRBBe  CONSTANTS 

• FOR  EACH  THERMISTOR 

COMPUTE  PLATINUM  PRPBe  TEHPE«AtU«E 
•COMPUTE  INDIVIDUAL  BRIDGE  RATtO 
COMPUTE  first  .'PPRBXImATIbsj  Tb  prT 
•CALL  CALLtNOAR 

COMPUTE  iTEPATIVrLY  A CORRECTION  TO  PRT  T9  YIELD  PINAL 

corrected  temperature 
Output  Data  for  validation 

,-RITF  BatCHNUM,DatE/PTPB8Be, veLTS,ZrR9RATl9, I N I TRAT I 0, BOOK, PAGE, 

N9MTEHP 

•fbr  each  thewh I STOR 

• WRITE  B*TwIDa  TmErmJd#rATI9#T68»UrES#RES 

•return 

HUELLER  BRIDGE 
•SUBRBUT I NE  HUELL 

•READ  TwEpHIsTbr  BATh  SEUuEnCF  NUMBERS  (ThErHID) 

«F0R  EACH  NOMINAL  TEMPERATURE 

•ReaDj  NBHINAL  TEmp  (N9HTEHP),DATE« BRtDGe  ZEPB  ERRSR  (ZER0ERP)j 
YbLtS  across  TmErh  I sT  9R  (VeLTS>*  s08*  anD  paGE  9F  qATA  j 
SCHLEICHER'S  NOTEBOOK,  INITIAL  BRIDGr  READING  WITH  N8RHAL  1 Ha 
CURRENT  THRBUGH  PRT  AT  MEaSUREmENT  TfMPERaTURe  (NRMLjMA),  AS  NRMLjMA 
BUt  R I Tw  CURRENT  WEVEHSED  <RVRs1ma).  SENSITIvITy  9F  INTFRP9L.ATING 
NULL  rRIDgE  RECORDER  IN  Bwhs/CHART  QtvISION  (COeFF),  pRBgE 

resistance  at  0 c *nd  0 M*  current  irzeRbi. 

•FbR  EACH  ThERMisTbR 

•rfaD;  change  in  Bridge  settings  fpbm  nphli^a  in  steps  of  .oooiohm 
(BRijnjha),  ahbunt  bridge  is  off  null  when  set  at  nrmljra^bri ji*a 
jm  ch^rt  divisions  (nullnima>,  ch4nGe  in  bridge  settings  frbh 

R VBS 1 h a IN  STEPS  9F  .0001  Ohm  (BRtjRiMAj,  AMOUNT  BRIDGE  IS  BFF 

null  whEn  set  at  rvrsiMa+bri jrima  in  chart  divisions  inullriHai, 
Change  IN  BRjDUe  SETTINGS  FR8M  RVpglAM*  IN  STEpS  9F  *0001  9HM 
( BRI jMi 4m* ) # AHBUNT  BRIDGE  IS  9FF  null  WHEN  set  at  RvRSiaHa* 

BR  I JR  i amA  IN  chart  DIVISI9NS  ( NULl.  Hj  A J , THeRH!ST8R  RESISTANCE 
(RES) 

• F9R  eaCh  THERMISTOR 

COMPUTE  PLATINUM  PROBE  T£MPERATURf  (TfeO) 

•COMPUTE  PRT  RESISTANCE 

• CflMPuTF  FirST  APPr9X  IMAT I9N  Tg  Prt 
•CALL  CALLENDAR 

COMPUTE  * CORRECTION  TO  FJRST  PpT  ITfRATIVELY 
butpUt  Data  for  validation 

•w«ITE  eATcHNuli,DATE,pTpR0RE,VOLTS»RZFRe,ZFR0ERR.C0FEE»B00<,PA3E 
NBMTEHP 

•FOR  EACH  THERMISTOR 

•write  bath  seuuENce  number. thermjo. prt  resistancEjptresi. TbA.RES 

•RETURN 

DIRECT  R.T  INPUT 
•SUBROUTINE  DIRECT 

•F BR  each  ThErhIsTbR 

• RF  AD  THERMID  AND  BATCHNUM 

•for  e*ch  nominal  temperature 

*rE‘D*  temperature  Amo  ThErhIsTbr  pFsf St*nCE 

•return 
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SECOND  UE  V£L  SUBROUTINES 

SUBROUTINE  C*LLFnDAR 

CaLL£NDAB  Cf>**PUTrS  A CoRRfCTION  TO  TWg  COMPUTED  PRT  TEMPEBaTUPE  3Y  AN 

iterative-  technique  specified  IN  THE  pmT  RRobe  hanual, 

•Compute  CBBBECn»N  (CORRl  BaSeO  ON  0EUta  <P<?T  CBNSTANt)  and  TE^EpATU«E 
lb  LESS  than  .00001 

•Compute  ANBTHfcH  CORRECTION  (DeLTAT)  BASED  ON  THUS  EAR  CORRECTED 

Te mper atuR*  and  sbhe  cbnstants. 

•FINAL  T68  waS  BBTH  CORRECTIONS  INCORPORATED, 

•RETURN 

SUBROUTINE  BFSCORR 

bescbbh  cbbbeCts  the  the9mistb«  pesistances<u«es)  eor  the  Resistance  of 
the  bath  thephistbw  leads  and  the  bridge  d j 4l  Resistances* 

•test  ebb  QHD.EU.aHEBBh  (E»  BRIDGE) 

tes-ob  t a lead  resistance  corrections 
NB.PRbCeeD  with  GUILDLJNe  BRJOGe  dial  CORRECTIONS 

• FOR  EACH  THERMISTOR 

»TtS7  EBB  URCS*Lt*100* 

TES-SET  CBBRECTIBN  EQUAL  TO  o 

nb-pbbCeed  with  cbhputing  correction 
*CA|_L  GETDIGITS 

^compute  correction  mr  each  dial  (Lead  , digits  onlyj 

• apply  TOTAL  CORRECTION  to  URfS,  OBTAINING  Res 

• For  E*Ch  TmFRHIsTbR 

•apply  lead  CORRECTION  to  RES 

•RETURN 

SUBROUTINE  GETDIGITS 

GETDIGITS  CONVERTS  a 6 DIGIT  'NUMBER*  MrO  A t DIMENSIONAL  ARRAY,  tSTORE' 
*£FR3  THC-  ARRAT  'STORE' 

*Onvekt  the  6 digits  bf  t he  resistance  Into  an  integer 
«PUT  the  DIGITS  of  the  INTEGER  RfSISTANCE  Into  'STORE* 

•return 

SUBROUTINE  ABCFIT 

ABCF I T Performs  a Lr  AST  SQUARES  fit  of  the  cal  I BRAT  I bn  data  TO  the 
function  v«a*h*x-*c*x**3 
-h£R£  Y» 1 /T68,  x»L0G(RES» 

AbCMT  is  RIJN  ONCE  MR  EACH  thermistor  IN  tmE  gATCH 
•f  PR  e*Ch  nominal  TEmpERatURE 
•COMPUTE  MEANS  of  X,  Y, X**3 

• For  each  NOmInA|_  temperature 

•subtract  me*ns  prbm  equation  to  be  MttfD  and  compute  summations 
necessary  for  DeteRmINjng  b,c.  Removing  Twr  means  Reduces  the 
roundoff  error. 

•Compute  BiC»*  CONSTANTS 
•for  EACH  nominal  TEMPERATURE 

•computf  estimated  values  bf  v and  t using  a,b,c  and  x,  and  residuals 
•return  ^ 

subroutine  Printout 

print  output  of  least  squares  fit  for  one  thermistor 

•PRINT  TWERMfO,  BatChnum 

•Print  Time  -of  C*lM  RUN< ICLk>'CALIBDATE*CA|  iBPi  ace 
•PRINT  A,B,C  CONSTANTS 

•RRjnt  y, yhat,ResIDY  For  eaCh  nominal  TEMpERATURr 

• RrJnT  rEs*T68.ThAT.i^EsIdT  F9r  Each  N0mInAl  TEmPFraTurE 
•return 
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C SUBROUTINE  PUNCHOUT 

C PUNCH  «,T  Data  PBINtS  and  ReSULtS  or  LEAST  SoUaReS  pit  FOR  each  theROJSTOR 
C RFBR  EACH  Batch 

c 8ATcHNgM,NynTHEHH«NuMNBHTEMP»CAL*B0ATE»C*Ll8PLACC 

c *cbhpute  nuhber  of  cards incardsi  required  tb  hbld  r,t  data  for  bne 

C THE«MISTB«  at  a PBINTS  Pe«  card 

c #FBR  EACH  THERHISTBR 

c •HU^CW  THERH!D«BATc**UH»CAltbDATE 

C #F BR  EACH  R.T  CARD 

c ^compute  the  subscript  bf  the  first  data  pair  bn  the  card  ipairii, 

C anO  the  Last  paiR(PairF) 

C # I F PAIKF  is  GREATER  TRAN  NUHNBHTEHp  then  set  PAIRF 

C EQUAL  tb  NURNBHTERP 

c rCbhpute  t«e  nuh»er  bf  r,t  data  pairs  tb  Be  punched 

C BN  This  card  (nbnCAroi 

C p PUNCH  NONCARD, NONCAR0  SFTS  bf  TftBfRES.  ANO  theRhid, BATCHNUH 

C RPUNCH  3 SETS  BF I 

C BF pR  EACH  ThERHJsTbR 

C *PUNCH  theRHID,BatCHNuh,CaL!BDatE#A,B#C#NUHNBHTEHR#RhSRESIDT 

C ^RETURN 


c PRBGRAH  CALIBB 

c 

C R.  PAYNE  9 JULV  1976 

c 


real  CALlBDATEO)#  CALIBPLAcEO) 

INTEGER  BATCHNUH,GHD»  ICL*C(*I 

call  tbdav(iclR) 

100  READ  ( 105#1000*END.150)  bATCHNUM,  qHD,  NUHThERH,  NUHNBHTEHp* 

• CALIBDATEiCALIBPLACE 

IF  (NUMTHERH,GT.25>  OUTPUT  .NUHTHERH  TB8  LARGE* I 
» Go  TO  900 

IF  { NUHN8HTEHP»GT • 101  OUTPUT  • NUHNBHTEHP  TBB  LARGE1! 

• oo  TB  9 oo 
IF  ( GMq*EQ« AnGyl L * C*ll  Oyl^o 
IF  (GHD.EQ.AHERBR)  call  guild 
IF  IGHd.EO.aHHueL)  call  hjelL 
|F  igmd.eq.ahoire*  call  oireCt 

DS  120  NTH.t#NUPTHERH 
CALL  ABCFIT 
call  printbut 

120  CONTINUE 

IF  CGH0.E0.AH0IREI  GO  TO  jqO 

call  punchout 

GO  TO  100 
150  CONTINUE 
900  CONTINUE 

1000  FORMAT  (T7»l3#  Tj6»A%»  T36*I?»  T5j,!g»  ?X#A?*  i*»*3»  jX*A g# 

• 1X,3a*» 


wsm 


SUBROUTINE  DIRECT 
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240 

2*1 

2*2 

2*3 

244 

245 
24fe 
2*7 

24% 

249 

250 

251 

252 
2^3 

254 

255 

256 

257 

258 

259 

260 
261 
252 

263 

264 

265 

266 

267 

268 
269 
2/0 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
2»2 
283 

2 44 

2«5 

286 

287 

288 
289 

'290 

291 

292 

293 

294 

295 

296 

297 

298 

299 


C 

••*••••*•• 

C SUBROUTINE  DIRECT  RpaDS  INpUT  TO  CaLIb  EROM  CapDS.  INPUT  IS  PROVIDED 
c as  values  er  o«ecise  temperature  and  resistance  eor  each  thermistor, 
c 

C Input;  CALIBB*n'IwThFrm,numn,jmTEmp 
C C.AW0S»TMERMI0,T68#RkS 

C OUTPUT;  CALIB«.TMERmIO# T68#«ES  PRO*  1 CALIBRATION  BATCH 
C 

C •.***4*.** 

c 

INTEGER  TMERMlD(Pb) 

DOUBLE  PRECISION  RES ( io* 25 ) * T6« ( lo* 25  ) 

00  200  NTM.^NUMTHERM 

Read  (105,2005)  THEH,1IO(NTH)  .BatChnUM 

READ  (105*2010)  (T68(NNT*ntm)#reS(NNT»ntwj,nnt.i#NUMN0mtemp) 

200  CONTINUE 
RETURN 
C 

2005  FORMAT  ( T j 3 » I 4 * Tp5» I 3 ) 

2010  P PRM a T (*(?P) > 

C 

* 

c 

SUBROUTINE  GUILD 

C 



c 

C GUILD  RtAOs  in  «*N  Cards  platinum  tweRMomeT^r  bpiDge  Data  and  THERMISTOR 
c Resistances,  correCts  por  leao  and  prioge  Resistance  (via  ReSCBRRj,  computes 
c TEmpeRatURes  prom  Ppt  Data  (via  CaLLeNDaR),  aNO  OUtpUtS  Both  input  and 
C Computed  Data  e*»r  validation, 

c 

C INPUT;  CALIBR.N'JMTWrRMjNUMNBMTEMPjBATCWNUM 

c LARDs-TRpR*in*PTPR*BEj^ERtjRATIft*NeMTEMP,  JN  j TP  aT  1 0*  VbLTS*  BoO<*  p A^E* 

c deltrat io*ures 

C OUTPUT;  CaLIHP.TMeRmi0,T68»RES(C0RReCteD) 

C PR INt.BatChnUM,Date*RtpRoBe<ZeRoRatIO, fNlTRATlO,VoLTS,B08K,pAGE, 

c N9MTtMP*T  HE  RMID, RATIO, T<,b#RES*UReS 

c 

c 

c 

REAL  N»MTEMp*DATE(3) 

INTEGER  B0e<*PA3E#PTPR0BE»BATHin(25  j 

double  p»fCision  ureS( io,25) , aLpha#delTa, initRatIO,Ze*oRatio, 

• RATI0(25)»DELTRATI0(25)*PT 

c 

C * i nput  * • • 

C 

READ  (105*21011 ( TH£R m ID(NTHj,nTH»j,num1 wE Rm ) 

READ  (105*2102)  PTPHOHE*  ZERORAT I 0 
DO  215  NNT" 1 * NUMNOMTEMP 

READ  (105*210*)  N(*mTEmP*OATE,  INITRATIO,vOLTS*bOOx*pAgE 
Read  ( i o5* 2106 ><batnid<ntm), deltrat io ( nth i,ures(nnt, nth), 

• NTW« 1 .NUMTWERM) 

C 

IE  (BATCwnUM.gT.51 1 CALL  RESCORP 


300 

301 

302 

303 

304 

305 

306 

307 
305 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 
326 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 
34  4 

345 

346 

347 
346 

349 

3b0 

351 

3*92 

353 

354 

355 

356 

357 

358 

359 


C 

C •••C»MPutE  TEmpfrATuRES  frOm  P^ATINuM  prOBf  tMPjjT 
C 

IE  (PTPROBE.EQ.#HBRAD)  ALPHA*. 00392577  j OELTA*j.l97 
IE  <PtPRbBe*eQ**hBU0Y ) ALpha*. 00 392636 4 1 DeLTA»1»4966 
IE  ( (PTPROBE. NE.4HBRAD)  .AND.  (PTPR0BE.NE.4Mfi'J0Y)  ) OUTPUT 

• • PTPROBE  INCORRECT* 

U8  21*  NTH*1,NU«THEHH 

Ratio (nth)* ini tratioydeLtRattoi nth) 

PT ■ ( (RATI9(NTH)/^eR0RATI0)*1»D0)/ALPwA 
CALL  CALLENOAR 

214  CONTINUE 
C 

C ***t)UTPUT*** 

C 

-RITE  (108.2190)  tjATCHNUH. DATE. oTpROoE. V0LTS. 

■ zeRoratio, initratio, book, page, nomtemp 

WRITE  (108.2192) 

-RITE (108* 2194) (BATHID(NTm,,theRmID,NTH,,RAT10(NTN), 

• T68(NNT,NTH),URtS(NNT,NTM),RES(NNT.NTH),NTH*i,NUMTHERH) 

215  CONTINUE 
RETURN 

C 

2101  EORMAT  (16H 

2102  EoRHAT  (T9,a4,  t25.E ) 

2104  EORM-T  (T^El*!.  TJ6.3A6.  T36.E.  T5?.E,  T64.lt  T75.I) 

2106  FORMAT  (TA,I2.  Tio*2E) 

2190  E ORH AT  < 1H1,9X, 'BatChnUM*',  j3/ 

• IOx.'OATE  DONE**.  3 A4/ 

• 10X. 'PTPROBE* ' «2X.  A4/ 

• 1 0*. ' VOLTS  ACROSS  TM£RM**,  E6.3/ 

• 10X, 'ZERORATIO*'#  E10.6/ 

• 10X. • INITRATIO**.  E10.6/ 

•lOX.'BOOX  NUMBER**.  15/10*. 'PAGE* *.  15// 

• 2X,*N0MINAL  TEmReR*TURe* ' . E6.2//) 

2192  EORMATj  * BATH  THERM* «Ti7. *P«T* .Tgg, i pRT* ,T35. *UNC0RRi.T*8. *CORR*/ 
a*  IDi. To, iNO.tTjfc. .RESIST*. T27, , TEMP,, T3*, .THERM  RES  THERmRES.) 

2194  EoRMAT  (IX,  *3,  J 7 , EU*6,  E9*4,  2ril*2) 

2195  EORMAT  (i  THERMISTOR  *.  U# * B*TCh  .«  I3,  ?X.3A4) 

2196  EORMAT  (4(E i2.4.E8.2») 

C 

* 

c 

SUBROUTINE  MUELL 

c 

* 

C SUBROUTINE  MUELL  READS  in  on  CARDS  pRT  bRIDgE  DATA  and  THERMISTOR 

c resistances,  computes  temperatures  erom  prt  data  <vta  callendarj,  and 
C 0U7PUTS  Both  ,NPUT  *ND  Comp UT£0  Daja  EoR  VaLiDaTiON.  ONL'  The  BRaDsHaw  PROBe 
C HAS  HEEn  U^eO  WITH  the  MUELLFR  bridge,  muell  Is  EOR  HISTORICAL  data  only. 

C the  MUELLER  BRIDGE  IS  NO  LONGER  USED. 

C INPUT!  CALI98,NUMTHeRM,NUMN8MTEmP 

c Caros-thermio/nomtemp,oate.zeroepr.voltr.boo<.paGE.nrmlima.rvrsi,1a* 

c RwRS14MA,C0EEE,RjER0.BRI JniMA,NULI NlMA.BRI JRlMA,NULL»lMA. 

c BQIJM14MA.NULLM14MA,RES 

c OUTPUT!  CALI8B*THERmID.T68.RES 

c PR J NT. B ATCHNUM, DATE. PTPROBE, VOLTS, RZERO.ZEROERR.COEEE, BOOK, PAGE, 

c N0MTEMP,THERMID,PT»ES,T68.RES 


HEAL  NR*Li"A,NULLNiMA(25),NULL*imA(»5>jNUU  m1*RA<25»* 

* ORJ Jn1*a(25> ,BRl JR1«A<25J,BR| JM1*Ma(25) 

DBUbLt  PRECISION  T8TALNiMA(25),reTAL'»lMA(2^,#MEANTBT,SHR# 

* PTHES<?5) 

•••inrut*** 

REAL)  ( 105*2*02>  <THERmID<nTH),nTW*i# NURTHERp) 

PTRRBBe^BBaO  1 aLPma* *00392577  j DeLTA*1*%97 
OB  2*5  NNT«1»NURN0NTERP 

READ  <105*2*061  '''8hTEmP*DATE»ZER8ERR*  V8|  TS«B88k*pAgE, 

% NRM(_1NA.HVRSlMA*RVHSl*MA,CBEF»r#RZER8 

READ  ( 105*2*08) (BhI jMlRA(NTH)#NULLNiMA(MTHj#pRI jRjmAjNTM)# 

$ NULLRiMA(NTH)#BHI JRi*RA<NTM),NULL*1*MA(NTW>,BES<NNT#NTH)# 

* NTW»1 »NUMTMERM> 

•••COMPUTE  TEMPERATURES  erbm  PPT  INPUT 
08  2AA  NTw,i,nUMThEWm 

TBTALN^MAtNTP) •NRMUlMA*BR I JN j MA ( NTW ) *NULLN jNA ( NTH ) *C8EFF 
TBTALRlMAl  ntR>*RvRS1mA*BRiOR1mA<NTPUnULLR1Ra<NtH>*CBeff 

meantbt.<T9TALNiMacNTW)*TBTALRiRA(NTw) )/2.Do 

SHR.Rv«SlAMA»RyRSlMA*BRI jP1AMA(NTW).pRI jRiMA(NTm 

* ♦(NULLMtAMA(NTH)/i.AD0«NULLRlMA(NTW) )*CBEPP 

PTHES^TH’-MEANTBT-SHR-’OOOlDO.iE^EP* 

NUN*  I NT (PTRES(NTH)  ) «20 

IE  ( (NUN.EO.S)tBR. (NUN.GE.7) ) PTRES(nTW).PTRES(NTN)«. 00001 
PT»(PTReS(NtR)-R2eb8,/(w^Ew8* ALPHA) 

CALL  CALLENDAR 
2**  CONTINUE 

•••OUTPUT*** 

-RITE  uo8»2A90>  BATCHNUMiOATE/pTpR9aE,v9LTS#RzER8, 

S ZeR8ERR#C9EPP*BBBK,PAGe#N9MTEmP 

-RITE  (108*  2*92 ) 

-RlTE  ( 108*2*9*)  (NTH*  THERM|D(NTw)»PTRES(NTw>/T68(NNT4NTm* 

» RES(NNT,NTH),NTH,i,NUMTMERM) 

2*5  C8NTINUE 

RETURN 

"2*02  P9RMAT  ( 161 ) 

2*06  E9RMAT  <tR,E,  t13*3a*,  t3*,E,  t52,E,  t63,I,  j7k,i/ 

* T9*Fi  T 26 . E / T**4E,  T59 * E*  T73.P) 

2*08  E9RMAT  ( T9.7E) 

2*90  P8RMAT  <1w1,9X#'BatChnUR»'*I3/ 

» 10x#  ' DATE  DONE ■ * * 3 A*/ 

* 10x* 'PTPR9BE* ' »2X,  A*/ 

*10X,'VBLTS  ACROSS  THeR"* * ,F* .3/ 

* iox* 'Resistance  be  probe  at  t.0  degrees. »*f9.6/ 

* 10**' ^ER8  WE  BRIDGE* «»E9. 6/ 

* 10x,'C9FrF.<BRMS/DlAL  Djv  at  1Ma’’'*f9.a/ 

» 10X«'B99<  NUMBER* Ii/lOX* 'PAGE  NUMbER*i*I5/ 

« '0'*9X* 'NOMINAL  TEMPERATURE* **E6,g//) 


56 


*2q 

*21 

*22 

*23 

*2* 

*25 

*?6 

*27 

*28 

*29 

*30 

*31 

*32 

*33 

*3* 

*3*5 

*36 

*37 

*33 

*39 

**0 
**1 
**2 
'**3 
*** 
**5 
**6 
**  7 


**8 

**9 

*50 

*51 

*•52 

*53 

*b* 

*5b 

*56 

*57 

*b8 

*b9 

*60 


*61 
*62 
*63 
*6* 
*6b 
*66 
*67 
*68 
, *69 
*70 


*73 

*7* 

*75 

*76 

*77 


*78 

*79 


2*«*2  Format  ('  BATH  TH£HM  PT  »ES  pt  TEmp  THeR*  RES') 

2*9*  FORMAT  (jx#  13*  17*  P11.6*  F9.*»  F11»2) 

2*95  P9R"AT  {•  THERMISTOR  •*  1*«  t patch  «,  I3,  2*«7A*) 

2*96  format  ( *(Ft 2.*,F8.2)> 

C 

••**•* 

c 

subroutine  CallfnDar 

c 

•••• 

C CAlLt^DA*  C9mpUTeS  A CORRECTION  t9  T Bv  aN  I TeRAT I Vr  TECHNIQUE  SpeC1pIeO 
C IN  the  HPT  P«9BE  HANUAL. 

c 

C input:  GUILD  99  HUELL-pT 
C 8UTPUTJ  GUILD  99  *UELL*T68 
C 

c 

c 

NC9UNT.1  , CORRl ■ 1 j C6P9.0 
TEHPaPT 

REpEAT  225,  PHJLE  OaBS(C9991.C99B)  ,Gt, 0.00901 

C9991.CP99 

t9P9«0ELTA»oi00*Ttpp*(  •OlOo*TEMP-l*Ool 

TEMP«Pt*C999 

NCBUNT.NC9UNT*! 

IP(NC*HNT.GT.i5)  OUTPUT  .CALUEN5A9  eUUAtIBN  OOpS  N9T  C9NVER3SU 
« RETURN 

225  CONTINUE 

DEL l*T»«o*50o,»9iOo#TEHP*  *»OlOo#TEMP*l»Op>*( tEHP/*i9*58DoM *00  * • 

« ( TfcHP/A30«7*00-l»00) 

T63(NNT,NTH>  aTEPP^DELTAT 

return 

C INPUT 



SUgROUT I Nr  9ESCPPR 
C 

* * 

C 9ESC93S  CORRECTS  THc  TM£RMISTOR  WeSISTANCES  r*®  TH£  RESISTANCE  8f  THE  gATM 
C LEADS  ANO  The  BRIOGf  DIAL  RESISTANCES. 

C 

C input:  GUILO  99  HUELL  • UH£S 
C g jtpu t 1 gujld  99  *ufll  • rEs 
c 

C 

c 

REAL  LtA09E5(2b)*BOCt0llO»lS3> 

iNTtQtR  st9re<<>»  *dig, dec  pals 

DATA  LEADRES/12( *3S)»12( «*0>/ 

0*TA  BOC/pOn*  iD-5*  20-5*  2<  0^0  >•?<  *5?*5  )*Pl*TV5)»  ?( -*0*5  )* 

« 2(0D0> 10-5, 2 ( 000 ) *ln«5,0D0,2D«5, 30-5,20-5, IQ. 5, 

* 6 ( ODO ) * 3( 1D»5  ) / 

C a-aER  B9I0GC  IS  USED  WITH  5q00  Sep!ES  TwEpPtST®RS.  NO  OlAL  CORRECTIONS 

c are  usrD  •* i tm  tmf  ep  hr i oge  put  le*d  Resistance  corrections  the 

c S a me  AS  p-OR  the  L and  N guarded  (.heaTSTpNE  g®  I OqE  • 

c 

IF  { GHD,Fn.*MERRRl  GO  TO  339 


^ * 
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*80 

*81 

*82 

*83 

*8* 

*85 

*86 

*87 

*88 

489 

*90 

*91 

*92 

*93 

*9* 

*95 

*96 

*97 

*98 

*99 

500 

501 

502 

503 
50* 

505 

506 

507 

508 

509 

510 

511 

512 

513 
51* 

515 

516 

517 

518 

519 

520 

521 

522 

523 
52* 

525 

526 

527 

528 

529 

530 

531 

532 

533 
S3* 

535 

536 

537 

538 

539 


C 

C • •«  CORRECT  EOR  B«IOOE  Of*L  RESISTANCES 
C 

OB  238  NTW«t,NU*THEH*i 
OCORR-O 

IE  (iJ»ES(NNT.NTH).LT.100»)  0C6R9«0  • j 
• 3*  T9  236 

CALL  GETOIGITS 
M(jLT«NU*DlGlT-2 
09  236  DIG'1.3 

DC999»0C8RR*B0C(STBRE(0lGJ»0fQ)*i0**MULT 
236  CONTINUE 

«ES(NNT*NTH) •uhes<nnt#ntn>*ocbrr 

238  CBNTINUE 
C 

*:  •••apply  le ao  resistance  corrections. 

c 

239  CONTINUE 

08  2*0  NTM«i ,NU*TMERM 

IE  (ReS(NNt,Nt«» .Lt.IOO.I  G8  TO  2*0 

IE  (QmD»EQ.*wERbR)  RES(NNT#NTMjbURES(NNt,NTM) 

RES(NNT,nth)«ReS(NNT#ntmj.leaDRes(NTM) 

2*0  CONTINUE 
return 
c 

C******* •#••#••••• 

c 

SUBROUTINE  OETOIGITS 

c 

•••••••••••••••••••«•••••••*••••••••••• ••• 

c 

c GETOIuITS  COVERTS  a 6 OloIT  * Number*  into  A j OImE^SIONAL  aRRAva  *STORe*. 

C the  values  BE  , STORE*  ArE  THe  DlOlrS  OF  »ES  STa^TINQ  with  the  highest  ORDER. 
C E0R  EXAMPLE.  R£S**17135  BeCOHeS* 

C STOWE ( 1 ) •*  STORE ( * ) «i 

C STORE * 2 ) » 1 STORE ( 5 ) «3 

C STORE*  3) «7  STORE*  5 ) "5 

C INPUTI  URES(NNT,NTH, 

C OUTPUTS  NUMDIGIT*ST8RE 

c 

c •••••••••• 

INTEGER  NUHniGlT#NUH#OIGlT*TRUNCNUH 
INTEGER  BLANK/*  '/ 

»EAi, 

c 

c •••zero  the  array  .store* 

00  21  NO  X • 1 « 6 

STORE (NDx) •BLANK 
21  CONTINUE 

R«URES*NNT,NTH! 

IE  ( R« LT • 1 00*  > r»0  TO  26 

C •••MAKE  twE  6 DIGITS  OE  the  RESISTANCE  INTO  AN  INTEGER 

0«R-100000 

oecihals-o 

REPEAT  ?3.  RH I LE  D.LT.O 
R*R*1 0 
0-R-100000 
JECIHALS«0ECIMALS*1 
23  CONTINUE 
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5*0 

5*1 

5*2 

5*3 

5**' 

5*5 

5*6 

5*7 

5*8 

5*9 

550 

551 

552 

553 
55* 

555 

556 

557 
55g 

559 

560 

561 

562 

563 
56* 

565 

566 

567 

568 

569 

570 

571 

572 

573 
57* 

575 

576 

577 


578 

579 

580 

581 

582 

583 
58* 

585 

586 

587 

588 
589 


590 

591 

592 

593 
59* 

595 

596 


597 

598 

599 


C Pur  INTEGER  INTO  STUPE 

NU**P 

DIGIT-7 

REPEAT  35.  WHILE  DIGIT. qT.o 
0IGJT*D1GIT»1 
TRUNCNUM»NUM/10 

STORE { OJrjI T )■  NUH.<  TRUNCNUM.iq) 

NUM»TRUNCNUM 

25  CONTINUE 
NUMO I G I T • 6 

26  CONTINUE 
RETURN 

C 

* • •• 

c 

SUBROUTINE  ABCFtT 

C 

* •*•••••••••«•••*•••••••**•.••••••••••••••••••••••••**•« 

c abceit  performs  a Least  squares  ejt  of  the  calibration  data  to  tme  function 
c Y»A*B*X*C****3 

c WHERE  V-l/T.  K«l&G(R> 

C 

C INPUTi  CALIBB.TA8.RES.NUMN0HTEHP 

C OUTPUT!  CaLI03.A,0,C, that, THAT. RES!Or,RfS!DT,BMSRESTDT 
C 

c *•• 

c 

REAL  RESIOV(io!.RESIDT(io).RMSReSIOT(25) 

double  PRECISION  T(10).XI10),TMeaN(xHeaN.X9HEaN<S0X0X,s0X0X3, 

• SOxDY.SOx30x3.SOx3DY.OY.Ox.Dx3. A( 25) <B( 35) «C( 25)  . YMATj 10) * 

• THAT(IO) 

C 

C ...COMPUTE  new  VARIABLES  AND  HeaNS, 

C 

YHEAN»XMEAN»X3MEAN»0. 

DO  302  NNT* 1 »NUHNOHTE*P 

Y(NNT),i.00/(T68(NNT. NTH) *273* 1300) 

X!NNT)»DLOG(RES«NNT.NTH) ) 

THE AN" ymEan* Y ( nnt  > /NUHN0HTEMP 
xMEAN*xHEAN*x ( NNT ) /NUHNBMTEHP 
X3riEAN.X3HEAN*(X«NNT)**3)/NUHN0*1TEHP 
302  continue 
c 

C ***SUBTRaCt  mEANS  TO  AVOID  ROUND  OFF  ERROR  !N  CURVE  fit 
C 

c •••COMPUTE  SUMS  FOR  CALCULATINQ  B AND  C CONSTANTS. 

c 

S0XUX«SOXOX3»SDX0Y«S0X30X3»SDX30Y.8, 

00  30*  NNT«! .NUMNOMTEMP 
Uv«Y(NNT)«TmEAN 
OX»X (NNT  j.XMEAN 
0X3"X{NNt)**7.X3MExN 
SDxDx«SD*Ox*Ox»Ox 
SDXDX3«SOXDXi*OX»0X3 
SO*Dr,SD*OY*oX#OY 
S0x30x3.S0x30x3*0x3*0x3 
S0X3OY.sOX30v*0X3»0Y 
30*  CONTINUE 

c 
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600 

601 

602 

603 

60* 

605 

606 

607 

608 

609 

610 
fell 
612 
613 
61* 

615 

616 
6 1 7 
618 

619 

620 
621 
622 
623 
62* 

625 

626 

627 

628 

629 

630 

631 

632 

633 
63* 

635 

636 

637 

638 

639 
6*0 
6*1 
6*2 
6*3 
6** 
6*5 
6*6 
6*7 

6*8 

6*9 

650 

651 

652 

653 

654 

655 


658 

659 


C ***C0MPuTE  A,  B,C  CONSTANTS 

C 

H<NTH)«(SDXOY«SnX3DX3-S0X0X3*SDX30v)/tSDXDx*SDX3DX3.SDXDX3**2) 
CfNTw)*(SOXOX«snx3ov.sDXOY*SDXDX3)/(SOXDX*sr)X30X3.SDXOX3**21 
A( NTH) .YMEAN.Bt  NTH) .x^EAN.Cl NTH)*X3MEAN 

c 

C •••COMPUTE!  pST  I MA  T£D  VALUES  OP  Y AND  T USK'3  A,B,C  AND  X 

C 

RMSRESIOT(NTH)«o« 

08  306  NNT«t*NUMNOMT£MP 

YMAT(NNT)»A(NTHUB(NTM)*x(NNTUC(NTH)*(xfNNT>**3) 

THAT ( NNT  1 ■ 1 .Dq/YHAT (NNT ) »p73 • 15Dq 
HESIDY(nnT1«YhaT<NNT)»Y(NNT) 

RES I OT ( NNT )»THAT(NNT)«T68(NNT#NTH) 

^^SReSJDT (NTH)«PMsHESIDT (NTH (♦RpSjOT (NNt )**2/NUMN0MTEMP 
306  CONTINUE 

pmshesidt(nth).sq»t (Rushes  jot (nth, j.iooo. 
return 


c 



c 

SUBROUTINE  printout 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


print  output  op  least  squares  pit  por  one  thermistor 


INPUTS  CaLI8B.theRm!D,BatCmnUM, ICLK.CaLIBDaTE# a»B,C. 

Y,  vhaT»REsIDY/REs<T68#THaT»PEsIDT.RmSPE'SI0T 
output:  print-same  as  input 


WR 

MR 

mR 

MR 

wR 

MR 

• 

mR 

mR 

WRITE 


ITE  (108**0001  ThERmID(NTH),bATChNUm 
ITE  (108**0021  ICLK*CALIBDATE*CALIBPLACE 
ITE  (108#*00*1  AtNTHl^BlNTHljCSNTM) 

ITE  (108**006I 
ITE  ( 1 08*  *008 1 

ITE  ( 1 08  , *0 10  1 ( NNT*  7 ( NNT  1 * VH  at  ( NNT  1 # RpS  I OV  ( NNT»_* 
NnT«i*NUMnOMTEMP) 

ITE  ( 1 08#  *012  1 

ITE  (108**01*1 (NNT,ReS(NNT,NTH),T68(NnT/NTH)#TWAT(NNT 
RESIDT(NNT>#NNT"l«NUMNeMTEMP> 

(108, *016)  RMSRESIDT(NTH) 


1* 


return 


*000  FORMAT  ( I iTHEPmistORi  I*  u*  I P*TCW I '»  !•*> 
*002  PORMAT  <'0  OATE  Of  COmPUTER  RUNS  *,*A*/ 

• '0  CALIBRATION  DATES'*  2x**2*  lx**3*  1X.*2// 

• ' CALIBRATION  SOURCE  I * * 3**1 

*oo*  pormat  t • o Calibration  constants’'' 

• 5x* • A • * * P15* 13/ 

% iO«**x,iS*'*  Pl5*13' 

« ' 0 ',**,' C*  ' * P15 .13 ) 

*006  PORMAT  ( ///25x* ' TAgLE  ®P  RESIDUALS') 

*008  PORMAT  ( »o' »T13» '0BSERvED'*T28* 'ESTIM 

• T i 6* ' Y ' # T32.'Y'*  T*8*'Y'l 
*010  eOrmAT  (IX*  15*  El5«10*  2e16,10) 

*012  PORMAT  (///Til* 'OBSERVED  OBSERVE^ 

• T 1 *, ' R ' , T26» ' T ' * T3  7* ' T ' * T*7*'T' 


ATED  ' » T*5* 'RESIDUAL'/ 
ESTIMATED  RESIDUAL'/ 


*014  format  nx<  15/  fiso/  pii**/  riots#  ru>ti 
*016  FORMAT  ( '0»*T33» 'RMS  RES  1 DUAL • , P7 .1, * MOEfl*) 

*022  FORMAT  (I*,  '/',l3z  '/'/At/  ""/At,  '•'/At.  '/'/Eil.Sz 

• ' / • /El  3 »9<  */•/ 12/  '/'/F5»2) 


'/'/Eu  .9/ 


SUBROUTINE  PUNCHOUT 


punch  Caros  "itm  r,t  calibration  data  for  piles  and  a,b,c  constants  for 

thermistor  INteRCOMPaRJSONS. 

I NPUT I Cal I BB { MUeLU . BATCHNUM, QMO, NU"TMfRM,  NUMN0MTfMP#  CAL I BOATEz  TMfRMl D, 

TftS/RES 

OUTPUT  I PUNCH.  SAME  AS  INPUT 


INTEOeR  NCaROS,PaIRI/PAIRP/NONCaRO 
IF  (QMO.EQ.*HERBR)  00  TO  50* 

WRITE  <106/5010)  BATCHNUMjNUMTHERM^NUMNOMTrMP^CALIBDATE* 

• calibplace 

NCAROS»NUMNOMTEMP/**l 
DO  502  NTH* i* NUMTHERM 

"RITE  <106,5025)  TMERMlO<NTH),BATCHNUM.«»ALtBDATE 
DO  501  NC«1/NCARD8 
PA  I RI ***NC«3 
PAIRP**»NC 

IF  (pAIRF.QT.NUMNOMTCMpj  pAIRF.NUMNOHTEMp 
N0NCaR0*PaIRF.PAIRI*1 

write  <106.5028)  NOnCaRO. <T6S< JP.NTWI .ReSI JP.NTH), 

• jP,PAIR1,PAIRF),TMERM1o<NTM),bATChNUM 

501  CONTINUE 

502  CONTINUE 
50*  00  5o5  «"l/3 

WRITE  <106,5035) «THERMlD<NTH>,BATCNNUM,eAL!BDATE#A<NTM># 

• b<nth),C(nth),numnomtemp/Rmsresidt<nth),nth*i,numtherm) 

505  CONTINUE 

return 

’5010  FORMAT  ( » BATCH# • , 1 3,  Tig/ 1 0M0*DIRECT i / Tf?. » NUMTHERM* ,, Ig, 

• T*0,'NUMNBNtfMP*',I2,  2X, a2,  IX, a3,  1*,a3/  3a*) 

5025  FORMAT  <*  THERMISTOR  i/I«,i  BATCH  */  13/  »x*A2z  iX/A3/  1X/A2) 
5o2»  FORMAT  ( N ( f 7* */  i X/ F5 *2/  iX)z  T71/I*/  i-t/TS) 

5035  Format  < I*,  ' / ' , 13/  */'/a2/  ' • » , A3,  **',a2,  '/',E13*»z  •/'/£!>•»/ 

• '/'/E13»R/  1/1/12/ 
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BaTCHNIUM*  67 
DATE  D6ME«  22  JUL  76 
PTPR©BE«  BU®T 
V«LTS  ACW3SS  THERM.  2*800 
ZER9PATI9-  2*555259 
INITHATI*.  2.555130 
B90K  NUMBER"  761 
PAGE*  *1 

NQMinaL  TEMPERATURE-  *00 


BATH 

THERM 

PRT 

P«T 

UNCBRH 

C®RR 

ID 

NO 

Resist 

temp 

therm  res 

therm  res 

1 

7 

2*555336 

•0076 

95046*20 

95045*15 

2 

32 

2*555337 

•0077 

?49?7 • oo 

949?5*  95 

3 

109 

2*555337 

•oo77 

95067*00 

95065*9? 

4 
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95799*10 

95798*05 

5 
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•OO77 
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94762*65 

6 
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•OO77 
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95350*05 

7 

163 
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•OO77 

94615*00 

94613*85 

8 

168 

2*555337 

•OO77 

94968*10 

94967*05 

9 

1*7 

2*555337 

•OO77 

95264  *6q 

95263*75 

10 

189 

2*555337 

•OO77 

95175*70 

95174*65 

11 

1*0 

2*555337 

•OO77 

95041 *5o 

95040*45 

1? 

23* 

2*555337 

•oo77 

95210*40 

95209*35 

13 

239 

2*555337 

• OO77 

97882*10 

9788i*2o 

1* 

2*0 

2*555338 

•OO78 

95042*40 

95041*30 

15 

242 

2*555338 

•OO78 

95495*00 

95493*9o 

16 

244 

2*555338 

•OO78 

95460*10 

95459*oo 

17 

319 

2*555338 

•OO78 

95115*00 

95113*90 

15 

320 

2*555338 

•OO78 

95188*40 

95187*30 

l9 

323 

2*555338 

•OO78 

95017*30 

95016*20 

2D 

324 

2*555338 

•OO78 

95167*30 

95166*20 

?1 

326 

2*555338 

• oo78 

95038*80 

95037*70 

2? 

328 

2*555338 

•OO78 

95012*00 

95010*90 

23 

345 

2*555338 

•OO78 

95l96*5o 

95l95*4o 

24 

357 

2*555338 

•OO78 

95189*00 

95187*90 

Input  data  verification  page  for  0°  C nominal  temperature 


TWcRmIST8«S  7 BATChI  67 
date,  be  compute”  runj  i4:o5  sep  2*/ *76 
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CALIBRATION  DATE!  23  JUl  76 

Calibration  sburcei  whoi 

Calibration  cbnstants 
A*  *000908784695* 

B»  *0002250005862 

C-  *00000011*9649 


TABLF  be  RESIDUALS 


BBSeR^ED 
Y 

1 *0036608908 

2 *0035952332 

3 *0035316941 

A *003A7o47lo 

5 *003*113163 

6 *0033541115 

7 *0032988069 


estimated  residual 

Y Y 

•0038A0«90l  -*0000000007 

•00359523**  *0000000012 

•0035316951  *0000000010 

•003A704679  -*0000000031 

•00'*M13182  *0000000019 

•0033541112  -*0000000003 

•0032988Q  7Q  *0000000000 


SbSERveO 

Observed 

ESTIMATED 

residual 

1 9 

R 

T 

T 

T 

1 

950*5 • 150 

• 007* 

*007B1 

•000051 

2 

74485*250 

4*9961 

4*99598 

•*ooooqo 

3 

58736*750 

10*000? 

10*0001^ 

* *0000R2 

4 

46648 • 750 

14*9953 

14*99559 

•000258 

5 

3728*. 250 

19.9920 

l9*99l8i 

* * 000163 

6 

29983*150 

24*99i5 

24*99i57 

•000028 

7 

24255*950 

29*9899 

29.98990 

••000002 

/ 

RMS  RESIDUAL  *126 

"DEg 


I ' HI  m I IHPPIff 


.......  ..  .... 
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Appendix  III 
CALCOLIN  Program 

The  CALCOLIN  program  is  used  to  evaluate  the  calibration  data  and 
compare  the  results  of  several  calibrations  on  a single  thermistor.  Its 
use  over  the  last  three  years  has  allowed  us  to  determine  drift  rates 
of  thermistors,  has  helped  us  to  detect  errors  in  calibration  data,  and 
to  evaluate  the  consistency  of  our  calibrations.  CALCOLIN  is  written 
in  Xerox  FORTRAN  IV. 

Following  the  listing  of  CALCOLIN  is  an  example  of  the  output  of 
the  program.  There  are  two  sets  of  numbers  for  each  thermistor. 

In  the  first  set  of  numbers,  the  first  four 
columns  are,  respectively,  thermistor  number,  calibration  batch  number, 
r.m.s.  calibration  residual  and  calibration  date.  The  numbers  in  the 
following  seven  columns  are  derived  in  the  following  way.  The  numbers 
at  the  tops  of  the  columns,  at  the  top  of  the  page,  are  temperatures 
( +0,  +5,  etc.)  in  °C.  The  numbers  directly  below  these  are  nominal 
resistances  for  this  thermistor  type  at  those  temperatures  from  one  of 
the  input  cards.  The  numbers  for  the  individual  thermistors  and  calib- 
rations in  Table  A for  each  thermistor  are  the  temperature  differences 
which  result  from:  (1)  calculating  a temperature  using  the  nominal 

resistance  at  the  head  of  the  column  and  the  A,  B,  C constants  for 
that  thermistor  and  calibration;  (2)  subtracting  the  temperature  at 
the  top  of  the  column;  (3)  multiplying  by  1000  to  get  the  temperature 
difference  in  m°  C.  The  absolute  value  of  one  of  these  temperature 
differences  has  little  meaning  except  to  tell  one  whether  an  individual 
thermistor  is  within  manufacturer's  specifications  or  not.  Comparing 
all  the  temperature  differences  at  a given  temperature  for  one  thermistor 
can  reveal  several  things.  If  the  numbers  stay  quite  constant  over  a 
period  of  several  years  as  in  thermistor  122,  we  can  say  that  the 
thermistor  is  stable  and  the  random  variations  in  the  temperature 
difference  give  us  a good  idea  of  the  repeatability  of  a given  tem- 
perature in  the  calibration  bath.  If  the  numbers  change  with  time 
as  in  thermistor  131  then  we  know  the  thermistor  is  drifting.  If  the 
temperature  differenc  * increases  with  time  this  indicates  that  the 
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thermistor  resistance  at  the  specified  temperature  is  increasing  and 
vice  versa. 

The  second  table  of  numbers,  for  each  thermistor, 
is  an  estimate  of  drift  rate  of  the  thermistor  together  with  data  which 
allow  us  to  evaluate  that  drift  rate.  Each  row  in  this  second  table  is 
the  result  of  computations  on  the  numbers  in  the  corresponding  column 
of  the  first  table. 

The  columns  represent: 

T - The  nominal  temperature  in  °C  as  in  the  previous  table 

N - The  number  of  calibrations  performed  on  this  thermistor 

TIME  - The  elapsed  time  in  days  between  the  first  and  latest  calib- 
rations of  this  thermistor 

YBAR  - Mean  of  the  temperature  differences  in  the  previous  table  for 
nominal  temperature  T over  all  the  calibrations 
SEMV  - The  standard  error  of  YBAR 
YINTER  -b  in  the  equation 
y = at  + b 

where  y = temperature  difference  from  Table  A 

t = time  in  days  from  the  first  calibration  to  the 

calibration  from  which  y was  obtained  and  a 
and  b are  the  result  of  the  least  squares  fit 
of  temperature  difference  vs.  time 
SLOPE  - a in  the  above  equation;  the  drift  rate  in  m°  C/day 

DRIFT  - The  drift  rate  in  m°  C/year;  SLOPE  * 365 

SELF  - Standard  error  of  the  least  squares  fit 

VSL  - Variance  due  to  the  slope  of  the  line 

VAL  - Variance  about  the  line 

CORR  - Correlation  between  temperature  difference  and  time 

SELF/SEMV  - An  indication  of  whether  the  least  squares  fit  with  time 
variation  is  better  than  a simple  mean 
CAL/VSL  - An  indication  of  how  important  the  variance  about  the 

line  is  compared  to  the  variance  due  to  the  slope  of  the  line 
- Following  the  seven  lines  at  nominal  calibration  temperatures 
is  a final  line  with  five  numbers.  These  are,  respectively, 
the  overall  means  of  YBAR,  SEMV  and  DRIFT,  the  square 
root  of  the  mean  of  VAL,  and  the  overall  mean  of  CORR 


The  final  table  is  a listing  of  all  the  input  cards. 
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Of  the  two  examples  shown,  thermistor  number  126  is  stable  with 
a standard  deviation  about  the  mean  of  ±1.2  m°  C.  Thermistor  number  76 
is  drifting  although  the  drift  rate  has  decreased  considerably  in  recent 
calibrations. 


I 


4. 

I 


1 

• 2 

3 

4 

• 5 
6 

7 

8 
9 

10 
1 1 
12 
13 
1* 

15 

16 
17 
1H 
19 
?0 
21 
22 
23 

’24 

26 

'27 

28 

29 

• 30 

31 

32 
r 33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
- 48 

49 

50 
■ 51 

52 

53 

• 54 

55 

56 

• 67 

58 

59 


7* 

c calcplim 

c calculate  delta  tevps  at  given  resistances  and  least  squares 

C PIT  HF  DELTA  TEMPS  T9  TI«e 

C* *Pp6 GR A mmER *R • p 4 Y^E * • 1 6 SEPTEMBER  1^76** 

DIMENSION  NTHER«(1000)ANH(  1000)  .DATE  (3,  1000)*  NPTSj  1000)* 
lpES  1 0 ( I'vrnl*  I ND  < 11  ),v,D(ipoo,»f'5(i0f)l-|)/DTE6Pf2«5<7)»DTEMP1  ( 35  ) # 

3n( ?5 ) , noav( 1 00 ) , I CL* <4)#oaYDIFE(25) 

QIU^LE  PRECISION  R£S(7)«,'RES47).NReSk(7)«T(7)» 

1 A < 1000) *2 < 1005) iC ( 1000 > a TEMP 

Data  T/373.i5#27R*i5,283»i5,288»i5,293»i5#2<48»i5#303,l5/ 

C NQA v THPMUSU  MATCH  a,9a  PAT  1 IS  1 jAN,  1973 

DATA  \DAY/3i/46*57#7?«66»o#o#0,0,OA254>?6g/a64»?67/g7g*?77*282» 

1 0J0,0'33l»352»37i»379»387»o<4oP»4n7A4l6#4S6*Qf477#0»5l5*0*528i 

? 5b7»564»n»673»7o2»7';3»73M*7o8#75o»754*795#8o8»8l6#822*0*9lS# 

3 927*946.973#9s»H*lo?6#o»llo9#  1 1 ?4#  oa  1 1 7}  , 0 , 1221 a q,  1 299*  1 33t a 

4 1349/ 

c******#***forhat  statements********** 

1010  format  ( T 1 0* i 1 > 

101?  FORMAT  (A4> 

im5  format  ( 14*  1 x/ 13#  1 XA  3 A3»  1 XA  3 ( Fj  3»9#  1 , I?j  1 y#F5,2«  9x#  2Ai> 

1020  FORMAT  ( 7 ( P 10*2) ) 

1110  F«RMat  (• 1 * ,5X, •THERMISTOR  CaLIpRatIBN  C0MP AR j 54* * , 4X, 4 A* > 

U12  FORMAT  (•  SLOpF  AND  DRIFT  AR£  MINUS  TH£  RATE  AT  WHICH  TEMPERATURE 
* W0JLD  AMPpAR'/*  TO  Change  IF  tw£  THERMISTOR  wrRf  mbuntfO  IN  A CO 
INSTANT  TEMPERATURE  hath') 

1115  FORMAT  (5ix» 'RESISTANCE/TEMRERATURE' ) 

1120  FORMAT  ( 4pX»  '4-0'  *<>XA  '45'  »5X#  '*1P»  #^X#  •♦i5*i5X»  '420'  »5X#  '425'A 
1 5X, '*30'  ) 

u?5  format  (4x#'Ther  days  bTCw  rsdl  Date  'a7F8.d 

1 145  FORMAT  ( • 1 ' ) 

1205  Format  < 'O'  ) 

1210  FORMAT  ( 1x»3A1»7I5»E6.2»1Xa3A3#7E8.i/5x*2A1) 

1?15  FORMAT  ( 1 9X j • T N TIME  THAR  SeMV  YINTER  SLOpE  DRI 

1 ft  self  vsl  val  C0RR  self/  val/») 

1216  pHRnAT  ( 1 ix  * ' c DAYS  mDEG  mOEQ  mDEg  mD/DA  mD/ 

1 yr  MDEG  M02  M02  SEmv  VRL'I 

1220  FORMAT  (?3X,F9.3»F7.3(20X#F9.3,  1 4X  , F 7 . 2 a F g, . 3 ) 

1314  FORMAT  { ' 1 ' » 4Xf  4A4 ) 

1315  FORMAT  ('  TH  HA  DATE  ABC 

1 PTS  RES' ) 

1320  fHRm*T  ( Ixa 3A 1 / 1 y » 2 1 4 » 1 x» 3A3a 3E 1 5 . 9# I 3* p5 • 2» 2x» 2Al ) 

£****«*• ***PROGRAM********** 

H a IfcND«0 
L • J*  1 

CALL  TBQAY( ICLK ) 

C**********  j nput  * 

C**  I Dn*0  USES  ONLY  G‘'sd  CALIBRATION'S 
C* • i dh* 1 uses  all  calibrations 
read  (105,1010)  I DB 

c.*  res  are  resistances  to  *e  substituted  to  get  temdpRatureS 

-f  READ  (105»1D?0>  (RES(  I ) a 1*1*7) 

C**wwjte  haGe  hea">ING 

"RITE  (IOSaUIO)  ICL< 

"RITE  (10*a1112) 

"RITE  «10*M115) 

"RITE  (108,1120) 

>" R I T E ( l 08a  1 12b  ) ( RES  ( I ) / I ■ 1»  7 ) 

« R I T E <10«a  1205) 

^no  continue 


*0 

*>1 

62 

63 

6* 

6b 

66 

67 

6* 

69 

70 

71 

72 
7 3 
7 4 
7 b 
7 6 

77 

78 

79 
30 

*1 

*2 

*3 

*4 

6b 

*6 

37 

<38 

69 

90 

•71 

9? 

93 

9 4 

9b 

9ft 

97 

93 

9q 

100 

101 

10? 

103 

104 
1 Ob 
106 
107 
103 

109 

110 
111 
112 
113 
11“ 

115 

116 
117 
116 
119 


c**  REAO  LARD  E»0M  CALlW  PRHGRAM 

c**«  0 IN  CrtLUM\  79  Mr  a NS  THE*"IST»«  WAS  BEEN  L9ST  *R  DrSTRBYEO 
c***  H IN  CHLIJWN  30  MrANS  HAD  CaLIBRatJBN 
c«  U 1 6 mENTS  1 E«W  EACH  NE  ft  a#b,c  CARO 

READ  ( 105. 16l5*e‘^D«405)  NTHERM(L),NB(L)»(DATE(I»L)#1»1#31# 
1 a<L>  ,6(L)  ,C(D/NPTS<L)  ,WESI0(L)  ^DIL)  ,MB(U 
IE  ( L»EQ» 1 > NTwrWMfo  ) .NTwe«M(  \ ) , NOAvi.NDAyiNgiL  ) ) 

I E { IDB»E0»  1 > G*  r»  ?O0 

I E (MO(L)  •Er3»lHO.0R»WB(U)»Ef3»lwP)  G*  TS  100 

C**********CALCUL  AT  IH NS ********** 

200  CONTINUE 

IF  (NThERM(L.).NE.NTHEWM(L-l)  ) G8  T»  4lQ 
210  CBN  T I NUE 

O*  CALCULATE  TEMPS  E»W  q I v E ni  res 
03  400  I » 1 / 7 

TLmp»1  ./(  A(L>*lML>*0Lf1G(WES<  I 1 > + C < L ) *0LBG  ( «ES  ( I ) )**3) 

0 TEMPI  j,  n,(TEMp-T(  I ) )*1000 

c**  0AY0IEE  is  N3.  me  Oats  erbm  EIRST  CALIB.  T8  CaLIB. 

c*«*  UNDfR  CmnSIOeRATIBN 

400  CONTINUE 

nA  yJ I EF ( J)  «NQA  Y ( NB ( C ) 1-NOAY1 

wRITE  (108/1210)  M0(L)#M9(L)«NTWERm(L),DAYDIEE(  J)»NP(LI# 
lRESIO(L)  , (0ATE<  I ,L) , I»l/3> /(0TEvp( J,M),H«1,7) ,M0(L>,m8<L> 
J»  J*1 
L-L  + l 
GO  T 8 100 

C*  J InCREmENTS  ! EBR  pACH  N£«  a,B,C  CARO  BUT  RESETS  when 
C***  T H£9m i stbR  NUMBER  CHANGES 
4Q5  I ENU ■ 1 
4 1 Q CONTINUE 

WRITE  (103/1215) 

-RITE  (103/1216) 

Jl« J*1 

C*#********LINEAR  EIT  °E  TEmp  DIEES**#**.***# 

vHaHM«DRTPTm« TtiARVARM«ALEM»CBRM«o 

OH  620  M«  1 , 7 

06  610  K ■ 1 * jl 
OTPMPj {<  )«0TEMp (K/M ) 

610  CBNTTNUE 

Tt.T(M)-?73*lG0 

CALL  LINETT  (nAYDIEE,DTE*pl/Jl/TE/v9ARVAR,VALE,C9R> 
YHaWM*  Y9  aPM  + Y«aR<7  7» 

ORIETM.OR  IETM4-0RIET/7* 

YbAHVARM.YBARVARM4YHlARVAR/(7.4Jl  ) 

ALEMp  aLPm*V ALE/7* 

C3WM.CMRM4CHR/7* 

630  CONTINUE 

Y -3  A HSU*. SORT  ( YHARVARM) 
aLSO"SORT ( ALEM ) 

-RITE  (108/1220)  YBAWM^YRAWSOfDRlETM, ALSO/CBRM 

-RITE  (108/1205) 

je  ( lENO'EQ* 1 ) G«  T 8 100" 

J»t 

\iOA  Y1  »NOAY  ( NB  ( L ) ) 

GO  T«  210 
1000  CONTINUE 

L»L”  1 

c**  «utuuti  imaGe  oe  input  cards 

-RITE  (103/1314>  ICL« 

•■RITE  (10«/1315) 


1?0. 

"RITE  ( 108, » 320 I (NO ( p ,«B( I >,NTWfcR*( I > ,N8( I > , 

121* 

1 ( DATE (K.t).K* 1.3). A(I),b(I>#C(I># 

122* 

? NRTsm»9Esiom*HDm#MB(i 

1?3. 

C 

1?*« 

c***« 

125* 

c 

126* 

SUBROUTINE  LINFIT  IX, Y,N,TEM»SDV2,VAL.C0RR> 

127. 

c 

128* 

c **•* 

129* 

c 

130* 

DIMENSION  X(25)»T(85),TIne(p5) 

131- 

c 

CONFUTE  *eanS 

132* 

XdARaYBAR.o 

133. 

09  100  I • 1 * N 

13%. 

XBAR«XBAR*X(  p/N 

135* 

YtJARa  vHAR*  y( I )/N 

136* 

1 DC 

CONTINUE 

137. 

C 

COMPUTE  SL^^E.^INTePCe^T^CORReLaTISN  (R), VARIANCE 

DUE  to  SLOPE  OF 

138. 

c 

LINE(VSL),  TOTAL  VARIANCE  OF  LINEAR  Fly(VLF) 

1 39* 

soxur»snx?*sDV2»o 

1*0* 

00  200  1*1 »N 

1*1« 

SJXUY»snx')Y*(»{I  ).XHAR)*(Y(  J ).YBaR) 

1“?. 

SUX2»S0X?*( X( t).XBaH)**2 

1*3. 

SUy2*SDy2*( y* I )• yOAR )**2 

1**. 

200 

CONTINUE 

1*5. 

slope«sdxdy/sdxp 

1*6. 

ORIFT»SLOPF*365 

1*7. 

*INT£R»Yf»AR-SLOPE*XBAR 

1*8. 

CeRN*SDX')Y/(  (SOX?*SOY2)**o.5) 

l*y. 

VSL"S0X0Y»*?/S0X2 

150* 

c COMPUTE  vARTANCf  Ab*UT  line  (VAL), STANDARD  E*ROR  OF  LINEAR  FIT  (SELF) 

1 5 1 • 

C 

AND  STANDARD  em°OR  OF  M^aN  VALUF  (SEmV) 

152. 

C 

IF  VSL  IS  MUCH  OREATFR  THAN  val  tnen  sloping  line 

IS  BETTER 

153* 

C 

PIT  TO  DATA  THAN  MEAN  VALUE 

15*  • 

VAL»0» 

155- 

DO  300  I ■ 1 » N 

156* 

vAL«V*L+{ y( I )-yINTER-SLOpE*X( t > )**2 

157* 

300 

CONTINUE 

158- 

SELF "SORT ( VAL'Nl 

159« 

SEMV»SU«T( SDY2/N) 

150* 

T 0 T I Mfc  »x(N)«xtl ) 

1 5 1 • 

RSLM*  y • 9 

162. 

IF  ( SUVR .Gy. #001  ) RSLM*SeLF/SE^V 

163. 

wvLM*0»o 

16*. 

IF  ( vSL‘GT.n.001  > RVLM-val/vsl 

165* 

£•••• 

166* 

-RITE  (108.3001)  TEM#N#T»TIME,ynAR,SEMv. VINTER, 

167« 

1 SLOPE. DR  I F T. SELF, VSL # VAL .CORK, RCLM,RVLM 

1 68* 

3001 

FORMAT  ( 11 X.F*.n, 13, I5#F9.3»F7.?#F10.3,F10.5#F9.3, 

169. 

1 F7. J,2F7.?,F7.*»2F8.3) 

170* 

RETURN 

171. 

END 

Tnfc^MjsTf1*  rALPWATIHM  ri?r-P4i?iS^M 
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